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DERIVATION OF LANGEVIN DYNAMICS IN A NONZERO BACKGROUND 

FLOW FIELD 

Matthew Dobson 3,4,1 , Frederic Legoll 2,3 , Tony Lelievre 3,4 and Gabriel 
(N ! Stoltz 3 ' 4 
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Abstract. We propose a derivation of a nonequilibrium Langevin dynamics for a large particle im- 
mersed in a background flow field. A single large particle is placed in an ideal gas heat bath composed 
of point particles that are distributed consistently with the background flow field and that interact 
with the large particle through elastic collisions. In the limit of small bath atom mass, the large par- 
ticle dynamics converges in law to a stochastic dynamics. This derivation follows the ideas of [3][6j[T] 
and provides extensions to handle the nonzero background flow. The derived nonequilibrium Langevin 

C^H, dynamics is similar to the dynamics in |19J . Some numerical experiments illustrate the use of the 

(— i obtained dynamic to simulate homogeneous liquid materials under shear flow. 
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*yv ' Molecular dynamics simulations have been increasingly used to bring atomistic accuracy to macroscopic fluid 
^^ ■ models [TUll!2tJll!2"Tll28"] . One example is the computation of the constitutive relation between the strain rate A = 
^— s ' Vu and the stress tensor <r(Vu, T) in complex fluids at temperature T, where one uses a microscopic simulation 
r»o . to determine the closure relation for the continuum equations |16] . We thus wish to simulate molecular systems 
t-H ' at temperature T that are subject to a steady, nonzero macroscopic flow, and the goal of this paper is the 
l~J derivation of a dynamics to sample such states. 
• ^h ' Simulation of molecular systems with a nonzero background flow is one goal of nonequilibrium molecular 
y\i ■ dynamics (NEMD) techniques [Hl§]. Several strategies have been proposed to sample molecular system under 
2 ' nonzero flow: for example, the SLLOD and g-SLLOD equations of motion [8,9,26,27 or dissipative particle 
dynamics |24j . Typically, these are used in conjunction with consistent boundary conditions such as the Lees- 
Edwards boundary conditions [S] in the case of shear flow or the Kraynik-Reinelt boundary conditions [2"5] 
for elongational flow. The basic equations of motion for these algorithms typically exhibit energy growth and 
need some form of modification if one wants temperature control for the system. Two common choices are 
the isokinetic Gaussian thermostat and the Nose- Hoover thermostat. It has been shown that the Nose-Hoover 
dynamics is non-ergodic for the NVT ensemble [171I18J , and furthermore we are not aware of an analysis of the 
ergodicity of the Gaussian thermostat for these types of nonequilibrium molecular systems. Instead, we work 
in a framework that leads to stochastic equations that have a Langevin-type thermostat. In the case of no 
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background flow, the standard Langevin dynamics has long been used to sample the NVT ensemble, and it can 
be shown to be ergodic. 

In the following, we derive a system of equations for a single particle in a microscopic heat bath with a 
nonzero, constant-in-time background flow in such a way that velocity gradient control and temperature control 
are incorporated at the same point in the model. The resulting system of equations, which we refer to as the 
nonequilibrium Langevin dynamics, is 

dQ = Vdt, 
MdV = - 7 (V - AQ)dt + adW, ^' 1 ' 

where (Q, V) are the particle's position and velocity, M is the mass of the particle, W is a standard Brownian 
motion, A is the homogeneous strain rate, and 7 and a are scalar constants satisfying the fluctuation-dissipation 
relation 

7=~<^ (1.2) 

where f3 = vp~ is the inverse temperature. This system of equations is a generalization of Langevin dynamics, 
which we recover in the limit A — > 0. In the sequel, we follow the derivation of Diirr, Goldstein, and Lebowitz 
(DGL) [6j[7] who consider the case of a heat bath that has zero background flow. 

We briefly summarize the ingredients of the mechanical model used in the derivation. A full description is 
given in Section 2. The microscopic mechanical model consists of a single large particle immersed in a bath 
composed of infinitely many small bath atoms. The mass of the large particle, M, is held constant while the 
mass of an individual bath atom, to, is a parameter of the mechanical model. The large particle moves in 
a ballistic trajectory until it collides with a bath atom, at which point an elastic collision occurs according 
to (|2.7[) below, causing a jump in the velocity. The heat bath is constructed in such a way that most of these 
jumps are independent events distributed according to a velocity measure centered around the background 
flow. More precisely, the bath atoms have a random initial velocity distribution that is centered around the 
desired background flow (|2.16|) - l|2.17|) . with its mean relative kinetic energy proportional to the macroscopic 
temperature. A microscopic dynamics for the bath atoms is chosen so that the velocity distribution of the bath 
atoms is preserved (one choice is given by (J2.18J1 ). up until they collide with the particle. A typical bath atom 
has velocity much larger than the velocity of the large particle, and such a bath atom will collide at most once 
with the large particle (a fact made rigorous in Appendix [Cj . The nonequilibrium Langevin dynamics is then 
derived as the limiting dynamics of the large particle when to — > (Thcorem l2.2p . 

The full description of the original DGL model as well as two approaches to incorporating background flow 
are given in Section [2] culminating in the main convergence result, Theorem 12.21 The proof of convergence of 
the heat bath model to the derived stochastic equations (II. ip is carried out in the Appendices. The proof is 
structured as in [7] , with added arguments to control how the flow in the heat bath affects the error growth. In 
Section [3] we include numerical results showing the application of the derived equations to computing the shear 
stress of a Lennard-Jones fluid. 

2. Models for non-uniform background flow 

In the following, we consider a system in R d (for d — 2 or 3) composed of a single, distinguished large particle 
immersed in a heat bath of light atoms that have mean velocity Aq at point q. We derive the equations of 
motion of the large particle in the limit as the mass of the individual bath atoms goes to zero. We note that 
we use the terms 'particle' and 'atoms' to differentiate the two types of objects, but the heat bath could well 
be composed of light molecules. We follow the arguments of Diirr, Goldstein, and Lebowitz (DGL) [7], who 
consider a single large particle placed in an infinite, constant-temperature heat bath with zero background flow 
(note that the case of a constant uniform background velocity is equivalent to the case of zero background flow 
by a change of coordinates). The only forces which act on the large particle are due to the heat bath (see 
Section [2.4.3l for extensions to multiple particles and more general interactions). In the limit of small bath atom 
mass for zero background flow, DGL recover the Langevin dynamics 

dQ = Vdt, 
MdV = --/Vdt + adW, ( 2 " 1 ' ) 
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Figure 1 . A large particle (of radius R = 2), surrounded by a heat bath with zero background 
flow, whose atoms are uniformly distributed in space. The velocity of an atom is distributed 
according to (|2.4[l , with /? = 1, m — 1, and A = yg. 

where W denotes a <i-dimensional Brownian motion, Q, V e R d denote the position and velocity of the large 
particle, and 7, a € R are scalar constants that are determined by the parameters of the large particle and 
the heat bath and that satisfy the fluctuation-dissipation relation (|1 .2[) . Both the drift and diffusion terms in 
the velocity equation are caused by the elastic collisions of the large particle with the bath atoms, and the 
randomness stems only from the random initial configuration of the bath. In this paper we extend the work of 
DGL [6l[7] to the case of a nonzero background flow. 

In Section r2.11 we review the construction of the heat bath and trajectory of the large particle in the zero-flow 
case of [7] before outlining in the remaining subsections two possible approaches for extending this to the case of 
nonzero background flow. The first approach, outlined in Section [2. 3[ applies to shear flows modeled by a heat 
bath with one or multiple unidirectional laminar flows. In this approach, the limiting nonequilibrium dynamics 
is not the dynamics given in (|1.1[) . For reasons outlined in Section 12.31 this first approach is not the one we 
focus on in this paper. We then describe a second approach in Section [2 . 41 involving a modified non-Hamiltonian 
dynamics for the bath atoms, which yields in the limit m — > the nonequilibrium Langevin dynamics (|1.1|) for 
general incompressible background flows. Note that, in contrast to the first approach, this second approach is 
not restricted to shear flows. 

In the following, we use bold notation for vectors and normal weight for scalars and matrices. Capital letters 
are used to distinguish the large particle's position and velocity (Q, V) from that of a generic bath atom (q, v). 

2.1. Heat bath with zero background flow 

Here we recall the model and results of 7 . The bath consists of infinitely many light atoms each with position 
q <G R d and velocity v € R d . All bath atoms have the same mass, m > 0, and zero radius. The initial bath 
configuration is drawn from a Poisson field (whose definition we recall below) with an ra-dependent measure 
given by 

d/im(v,q) = A m / m (v)dqdv, q,veR d , (2.2) 

where 

A m = m-^ 2 X 
is the expected number of atoms per unit volume and f m is the scaled probability distribution on the velocities. 
The velocity probability distribution scales like 

/ m (v)=m d / 2 /(™ 1/2 v), 

which means that for a single bath atom, the expected speed E m (|v|) = J Rd |v|/ m (v)dv is proportional to 
m -1 ' 2 . This scaling ensures that the average kinetic energy per atom is constant in the limit m — > 0. The 
distribution /(v) is assumed to be rotationally invariant. The quantities 

$i = l I |«i|*/(v)dv, for* = l,...,4, (2.3) 

^ .lad 



where v% denotes the first component of v, are assumed to be finite. For a set S C M. d x R d in phase space, 
with measure denoted by fi m (S), the number of atoms in S, denoted N m (S), is a Poisson random variable with 
parameter /i m (S), so that for k > 0, 

F(N m (S) = k) = ^pH e -^(s)^ 
k\ 

Figure Q] shows one possible realization of the heat bath with the choice 

f(v) = Z- 1 exp(-^\v\ 2 Y (2.4) 

(•> \ d/2 
where Z = I =£■ 1 is the normalization constant. 

As soon as the initial condition of the heat bath process has been chosen, the evolution of the system is 
deterministic. The bath is initialized and a large particle with a finite radius R and mass M is placed in the 
bath at (Q(0), V(0)), where we note that the initial condition is independent of m. Any bath atoms that are 
underneath the large particle at the initial time, that is, those such that |q(0) — Q(0)| < R, are removed from 
the bath (we recall that the bath atoms have zero radius). The bath atoms have no self- interaction and, aside 
from collisions with the large particle, obey the dynamics 

da = vdt, 

rf v = 0. ™ 

Likewise, the large particle's position and velocity evolve according to 

dQ m = V m dt, 

dV m = 0. ^ 

We explicitly denote the m-dependence of the mechanical process, while for notational convenience, we omit 
the m-dependence of the bath coordinates. 

When a bath atom encounters the large particle, it undergoes an elastic collision. For a single collision, let 
e n denote the unit vector in the direction from the bath atom to the particle center. Let V n = (V • e n )e n and 
V t = V — V n denote the velocity of the large particle in the normal and tangential directions. Similarly, define 
v n = (v • e n )e n and v t = v — v n for the bath atoms. We distinguish between the component of V in the normal 
direction V n = (V • e n )e n and its magnitude V n — V ■ e n . The elastic collision rule is 

V{ = V t) v( = v t , 



M — m , , 2m 

V n H 

M + m M + i 

M -m 2M 



V ' n ~ MT^ Vn+ TTT^ Vn ' (2-7) 

-V n , 



M + m M- 

where primes denote the after-collision velocities. The generated trajectory (Q m , V m ) is called the mechanical 
process. It is possible to only consider the case of single collisions since it has been shown that pathologies 
such as multiple simultaneous collisions or an infinite number of collisions in a finite time interval have zero 
probability [71 Appendix] (see also Appendix [EJ. As a result of the above, the mechanical process is well-defined 
on the finite interval [0,T]. 

Each initial condition of the heat bath corresponds to a realization of the mechanical process (Q m , V m ). 
The velocity V m is defined to have right-continuous jumps so that the trajectory is a cadlag function, that 
is, a piecewise continuous function that is everywhere right-continuous and has a left-hand limit. For a fixed 
initial condition of the heat bath, the particle trajectory (Q m , V m ) is a deterministic process; however, since 
the bath's initial condition is a random variable, the particle's position and velocity are also random variables. 

Remark 2.1. The trajectory (Q m , V m ) is not a Markov process. It would only be Markov if the rate of collisions 
with bath atoms did not depend on the history of the particle; however, this is not the case because of two 



effects. First, recollisions are possible with bath atoms that are moving sufficiently slowly. Second, certain 
collisions, called 'virtual collisions', with slow- moving atoms are made impossible due to the particle 'sweeping 
out' a path behind it. The large particle creates a wake behind itself, and whenever it decelerates, there is 
a minimum delay before slow moving atoms can hit it from behind. To show convergence of the mechanical 
process to the limiting Langevin dynamics (I2.1[) . one must show that these history effects are negligible in the 
limit m — > 0. 

We now define the appropriate topology in order to precisely state the convergence result of [TJ. We fix 
a finite time T and let T>([0,T}) denote the space of cadlag functions on the interval [0,T]. For functions 
f,g E T>([0,T]), we define the Skorokhod metric (see for example [2j) 

0"sk(/,sO = inf max{||A-Id||ioo( [ o,T]),||/-S , °A||ioc([o,T];R<i)} ) (2-8) 

where Id : [0, T] — > [0, T] denotes the identity map and A is the set of all strictly increasing, continuous bijections 
of [0,T]. For vector- valued /, we apply the Euclidean norm in space so that the L°°-norm above is given by 

ll/IU°°([0,T];R<*) =eSSSUp|/(»|. 

te[o,T] 

The mechanical process (Q m , V m ) defined above is in 2?([0,T]). 

We now recall the following definitions for the convergence of random variables. 



Definition 2.1. Let Z m and Z be random variables in the metric space (2?([0, T]), <7 s k). We denote convergence 
in law by (Z m )r 0<t<T -i — ^{o<t<T}> that is, E/(Z m ) — > E/(Z) asm->0 for all bounded, continuous functions 
f : T>([0,T]) — > R. We denote convergence in probability by (Z m )r 0<t<T j — K-Z){o<t<T}; that is, for all e > 0, 



lim F({a sk (Z m , Z) > e}) = 0. 

m— fO 

The following theorem gives not only convergence in law of the mechanical process to the Langevin dynam- 
ics (|2.1|) . but also an explicit expression for the coefficients of the limiting dynamics: 



4Ai? d - 1 S d _ 1 



*i, 



; $3 



(2.9) 



where Sd-i is the surface area of the (d — l)-sphere E> d 1 , $i is defined in (|2.3p . R and M are the radius and 
mass of the large particle, and A is a density parameter of the bath in 



Theorem 2.1 ( [7J Theorem 2.1]). Consider the mechanical process (Q m , V m ), defined in (|2.5[) , (|2.6p . and (|2.7|) 
with bath measure defined in (|2.2p . and let (Q,V) denote the solution to the Langevin equation (|2.ip with 
coefficients (12.91) . For any T > 0, 

r 

(Qm, V m ){ <i<T} KQ) V ){0<i<T} 

as m — > 0, where convergence is with respect to the Skorokhod metric (|2.8[) on 2?([0, T]). 

The resulting dynamics (|2.ip satisfies the fluctuation-dissipation relation (|1.2p provided that 

$ 3 2' 

which holds true for /(v) as in (I2.4p . for example. In this case, the temperature of the large particle is identical 
to that of the bath. The main result of this paper, Theorem 12. 2\ is the generalization of Theorem 12.11 to the 
case where the heat bath has a nonzero background flow. 



2.2. Generalization to a nonzero background flow 

We generalize the model described above by considering a bath measure d/i m (q, v) where, in contrast to 
the previous section, the initial distribution of the velocities depends on the position. We give two possible 
approaches to the generalization. In the first, which applies only to shear flow, the velocities are restricted 
to laminar profiles. For the second, which applies to general incompressible flows, we change the microscopic 
dynamics that the bath atoms satisfy in order to preserve the distribution of atoms. In both cases, we require 
that the bath measure and dynamics satisfy the following hypotheses: 

The average velocity at point q* is equal to Aq* in the sense that (HI) 

/ / v/i m (dq,dv) 

for all q* e K d , lim — j — q ,£ exists and is equal to ^4q*, 

^ / / fi m (dq,dv) 

JR d JB(q*,e) 

The measure d/i m (q, v) is invariant under the bath dynamics in the absence of collisions, (H2) 

where 23 (q, e) denotes the ball of radius e about q. The invariance of the bath measure means that, except 
for the effects described in Remark |2.1[ the large particle experiences collisions with a time- independent rate. 
We note that satisfying (|H1|) and (|H2[) is not automatic. For example, if the velocities v have a Gaussian 
distribution centered around Ac\ and the atoms follow ballistic trajectories (|2.5[) , then the system satisfies (|H1|) 
but not (|H2|) . In both generalizations described below, the bath atoms interact with the large particle only 
through elastic collisions as in (|2.7p . 

In Section l2~3l we describe a laminar flow model for the specific case where the background flow is a shear. All 
bath atoms are restricted to move in one of the coordinate directions, following otherwise ballistic trajectories. 
The resulting dynamics is not the same as the nonequilibrium Langevin dynamics (jl.l[) that form the main 
focus of this work. In Section 12.41 we do not restrict the velocity directions, and we choose a modified (non- 
Hamiltonian) dynamics for the bath atoms which leaves the chosen bath measure invariant for any traceless A. 
This approach leads to the nonequilibrium Langevin dynamics (jl.ll) as limiting equations for the large particle 
in the limit m — > 0. 

2.3. Laminar flow models 

For this subsection, we restrict ourselves to a heat bath in M 3 under shear flow with the specific strain rate 



A 



s 





(2.10) 



for some given s G R. We describe two variants of a laminar flow model, which both create a background shear 
flow by restricting the atom velocities to be parallel to the coordinate axes. 

2.3.1. Single laminar flow 

We enforce a steady shear flow for the bath atoms by choosing initial velocities that are nonzero only in the 
first coordinate, ei = [1,0, 0] T . We choose the initial configuration of the bath atoms according to a Poisson 
field with measure 



d/i m (q,v) = X m Z 1 m 1/2 exp f --m(v 1 - sq 2 ) 2 j S(v 2 )5(v 3 ) dqdv, 



(2.11) 



1/2 



where Z = ( =£■ J and S denotes the Dirac distribution (see Figure [2] for an example initial condition) . The 

bath atoms follow ballistic motion (|2.5j) . and it is easily verified that the bath satisfies (|H1[) and (|H2|) . The bath 
atoms undergo elastic collisions (|2.7j) with the large particle. 



Figure 2. A large particle (of radius R = 2) surrounded by a laminar heat bath whose atoms 
are distributed according to (|2.11| . The velocity is centered around the constant shear profile 
Aq, where A is given in (I2.10[) . Since the bath atoms may only flow parallel to the shear, the 
variance around the mean velocity along each flow line 52 = c is constant for all time. We have 
chosen s = 0.1 /3 = 1, m = 1, and A = j^. 



Ignoring the effect of recollisions, we formally compute the limiting dynamics as m — ¥ to be 



dQ = Vdt, 
MdV = - 7 (V - AQ)dt + adW, 



(2.12) 



where 




lV2^XR 2 I 3 



lp 
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The details of this computation are given in Appendix |FJ We note that for the laminar flow model considered 
here we have not rigorously proven a convergence result and, in particular, we have not justified ignoring 
recollisions. On the other hand, a proof is carried out for the model presented in Section [2^1 

The resulting dissipation and noise terms are both anisotropic, which is due to the fact that the background 
laminar flow is itself anisotropic. These terms are larger in the direction of the flow, and the dynamics do not 
satisfy a standard fluctuation-dissipation relation in which ^^ 1 aa T = 1, since 




(2.13) 



Furthermore, in the case where the shear flow is identically zero, A = 0, the dynamics does not reduce to the 
standard Langevin dynamics (|2.1[) which was derived for the case of zero background flow. Since we desire a 
system of equations that reduces to Langevin dynamics in the case A — (in particular because the Langevin 
dynamics is a widely accepted dynamics to model fluid flows at equilibrium), we consider a modification for 
removing the anisotropy. 



2.3.2. Multiple laminar flows 

One may attempt to create an isotropic flow by overlaying multiple background flows. We consider the 
superposition of three distinct laminar flows, each one with all velocities aligned in a single direction, e^. Then 
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Figure 3 . Two components of the laminar bath which are superimposed to create an envi- 
ronment for the particle consistent with the measure (|2.14j) . We have chosen s = 0.1, /3 = 1, 
m = 1, and A 



16' 



we choose the initial bath coordinates according to a Poisson field with distribution function 



d/i m (q,v) 



-A m Z 1 m 1/2 exp f -—m(v\ - s<? 2 ) 2 ] <5(«2)<5(v3), 
+ -\ m Z~ 1 m 1/2 5(v 1 ) exp f -~mi| j 5(u 3 ) 

exp I --^mv 3 



(2.14) 



+ iA m Z- 1 m 1 /^(z; 1 ) ( 5( U2 ) '' '' - 



dqdv. 



Each atom moves in one of the three coordinate directions, since the presence of the delta distributions means 
that at most one of the velocity coordinates is nonzero. Since the bath atoms do not have any self-interaction, 
the different flows do not 'mix' in any way, and in particular this bath measure is invariant under ballistic 
motion. Figure [3] depicts a realization of this bath in the two-dimensional case. 

According to the formal calculations carried out in Appendix|Fl where we again ignore the effect of recollisions, 
we find that the limiting dynamics of the large particle is 



dQ = Vdt, 



1 



where 



4y/2^XR 2 



AQ)dt + ddW, 



-1 1/2 



(2.15) 



7 = 



VP 



4y/2^XR 2 



ip 



As desired, the dynamics in (|2.15[) is no longer an anisotropic equation in contrast to (|2.12[) since 7 and a are 
scalar quantities. We note that the damping term on the large particle is relative to half the velocity of the 
background shear flow and that the fluctuation-dissipation relation is not satisfied, since |7 _1 u 2 = | 5^ 1. While 
we have chosen the inverse temperature for each bath equal to j3 in equation (|2.14[) . we find in Appendix IF1 that 
there is no choice of inverse temperatures for the three laminar bath components that allows us to construct an 
isotropic damping in which the damping term is relative to the full background flow. We note that this issue is 
not specific to the three dimensional case considered here. 

We do not pursue this model further for two reasons. First and foremost, it is not clear how to generalize 
these laminar bath flow models from simple shear flows to more general incompressible flows. Second, a heat 
bath with temperature proportional to /3 _1 and strain rate A gives a macroscopic dynamics whose temperature 
is proportional to \ji~ x and whose damping is relative to \A. We expect damping relative to A as in (|1.1[) , 
which agrees with NEMD dynamics such as g-SLLOD (see Section 12.4.21 for a description of the g-SLLOD 
dynamics and their relation to (jl.lj) ). While it would be possible to choose bath parameters Abath = 2 A and 




Figure 4. A large particle (of radius R — 2) immersed in a heat bath whose atoms are 
distributed according to (J2.17I) . The velocities have Gaussian distribution centered around a 
shear flow, with A given in (I2.10[) . We have chosen s = O.f, /3 = 1, m = 1, and A = i. 

/3bath = \$ to give the desired macroscopic parameters, A and /3, to the large particle, we instead consider 
a different approach which both handles general incompressible flows and shows agreement between the bath 
parameters and the parameters of the derived dynamics. 

2.4. Non-Hamiltonian bath dynamics 

In this section we describe an approach based on modifying the bath dynamics to no longer follow bal- 
listic trajectories (|2.5[) . The particle velocity accelerates so that the distribution of velocities relative to the 
background flow is preserved. 

2.4.f. Model and convergence results 

We consider a bath whose initial condition is characterized by the measure 

d^ T „(q, v) = A m / m (q, v) dqdv for q, v £ R d , (2.16) 

where f m has the form 

/ m (q,v)=m d / 2 /(™ 1/2 (v-Aq)) (2.17) 

for A £ R dxd (Figure 0} and where A m = m _1//2 A. As in Section I2TT1 we assume that /(v) is a rotationally 
invariant probability distribution, and that the first four moments are finite. Any distribution of the form (|2.I6[) 
with (|2.I7|) satisfies (jHIj) . We additionally assume that / is a decreasing function of |v|. This additional 
assumption is used below when approximating the trajectory of the mechanical process with a Markov process 
in Appendix ID1 

In order to satisfy (|H2|1 . we now change the underlying bath dynamics and consider the following non- 
Hamiltonian dynamics for the bath atoms: 

rfq = vdt, 

dv = Avdt. (2 - 18) 

It is still the case that bath atoms do not interact with one another, but now the bath atoms accelerate in a 
way that leaves the desired velocity distribution (|2.16j) - (|2.17j) invariant. The corresponding Liouville equation 
is 

dtp(q,v,i) = -V q • (\p(q,v,t)) - V v • (Avp(q,v,t)). 

As long as tr A = 0, which is equivalent to having an incompressible background flow, any function of the 
form ,o(q, v,i) — f(v — Aq) is invariant under the dynamics. This invariance is our motivation for choosing 
the dynamics (J2.I8I) . and in particular, we are not restricted to shear flows. We interpret the dynamics by 
considering the relative velocity 

v = v-,4q. (2.19) 
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In terms of (q, v), the dynamics (|2.18p is 

da = (Aa + v)dt, 

d!Lo. (2 - 20) 

Thus the relative velocity, v, does not change in time and the velocity of the atoms is at all times equal to 
the background flow plus an initial perturbation of thermal origin. However, the choice of dynamics (|2.18j) is 
motivated by the fact that it satisfies (|H2[) . rather than by the above physical interpretation. 

As before, the large particle undergoes elastic collisions with the bath atoms according to the rule (|2.7[) . In 
the limit m — > 0, the particle dynamics converges to the nonequilibrium Langevin dynamics (|1.1|) which we 
recall here: 



dQ = Vdt, 
MdV = -7(V - AQ)dt + <rdW, 



(2.21) 



where, as in the zero background flow case 

4A£ d -%_i, r4Ai? d - 1 5, ; , 

7 = ~ d *i, 



§ 



(2.22) 



where Sd—i denotes the surface area of the (d — l)-sphere, and $i is defined in (j2.3p . As in the case of zero 
background flow in Section 12.11 this dynamics satisfies a standard fluctuation-dissipation relation (jl.2p , with 
temperature equal to the bath's temperature, provided that 

$ 3 2' 

which holds true for the Gaussian distribution /(v) = Z^ 1 exp(— lr|v| 2 ), for example. We establish the following 
result in the Appendices IATIE1 

Theorem 2.2. Consider the mechanical process (Q m , V m ), defined in (|2.6p . (J2.7I) . and (J2.18P with bath measure 
defined in (12.161) and (|2.17p . and /et (Q, V) denote the solution to the nonequilibrium Langevin equations (|2.2ip 
wii/i coefficients (|2.22[) . for any T > 0, 

(Qm, V m ){ <t<T} KQ) V){o<*<T} 

as m — > 0, where convergence is with respect to the Skorokhod metric (|2.8p on X>([0,T]). 

The structure of the proof is the same as in [7] : we use a Markov approximation to the mechanical process and 
split the convergence proof into two steps. In Appendix [2] we define the Markov process that approximates the 
mechanical process and outline in more detail the convergence proof. In Appendix [B] we show that the Markov 
process converges in law to the solution of the nonequilibrium Langevin dynamics (|2.2ip . In Appendices [Q and [Dl 
we show that the difference between the Markov process and the mechanical process converges in probability to 
zero. These two convergence results are enough to deduce that the mechanical process converges in law to the 
nonequilibrium Langevin dynamics [5] . In Appendix [El we show that the mechanical process is well-defined for 
almost every initial condition up to a positive stopping time. 

2.4.2. Relationship to g-SLLOD 

The g-SLLOD equations of motion [51127] (also called p-SLLOD [S]) are given by 

dQ 1 = (AQ. t + V,)dt, 

(2 23) 

MdV i = (-V Qi E(Q)-MAV i -MAAQ i )dt. 

Note that these equations are written in terms of relative velocity, Vj = V$ — AQ^ . This system of equations is 
used to simulate a molecular system in a non-zero background flow. The particles interact through the potential 
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field -E(Q). The g-SLLOD equations are simply a change of variables from the Newton equations of motion in 
the reference frame 

dQi = V, dt, 
MdV % = -V Qz E(Q)dt, 

to the relative velocity. If we add a Langevin damping term and noise term to the dynamics in (|2.23p , we have 

dQ i = (V i + AQ i )dt, 
MdVi = (-V Qi J5(Q) - MAVi - MAAQ t )dt - 7 V, ; dt + adW. 

Transforming back to the reference frame, this gives 

dQ t = Vidt, 
MdVi = -V Qi E(Q)dt - 7 (V 4 - AQ t )dt + adW, ^ '" M > 

which is (|2.21[) . with the addition of an external potential E. Thus, our construction is consistent with the 
application of a Langevin thermostat to the g-SLLOD equations of motion and thus provides a derivation of 
the g-SLLOD equations from a heat bath model. 

2.4.3. Generalizations 

We have chosen the microscopic dynamics (|2.6|) for the large particle and (|2.18p for the bath atoms. First, 
we may make a more general choice of microscopic dynamics for the large particle between collisions, 

G'^lrii * m Ob, 

Mdv m = F(Q m ,y m )dt, f2J;>! 

replacing (|2.6[) with (|2.25p in the definition of the mechanical process. We assume F is uniformly Lipschitz 
continuous on M. d x M. d and Vy • F(Q, V) = which implies that both the mechanical process for finite m and 
the limiting SDE are well-posed (the proof given in Appendix lEl for the F = case depends on the fact that the 
flow map of the mechanical process preserves Lebesgue measure). The limiting dynamics of the large particle 
becomes 

rfQ = Vdt, 
MdV = F(Q, V)dt - 7(V - AQ)dt + adW. 

The proof of this result is a straightforward extension of the analysis given here, where the Lipschitz continuity 
allows us to bound the particle trajectory over finite time intervals. 

For the specific case of F(Q,V) = MAY with tr^4 = 0, which is the same acceleration term as in the 
microscopic dynamics chosen for the bath, we find the stochastic dynamics 

dQ = Vdt, 

MdV = MAVdt - 7 (V - AQ)dt + adW. ^' 26 > 

This dynamics was obtained in |19j for the case of shear flow by applying the Kac-Zwanzig formalism, thereby 
obtaining a generalized Langevin dynamics, and assuming that the memory kernel converges to a delta function 
(that is, there is no memory in the system). The dynamics (|2.26p is of particular interest since we can explicitly 

check that p(Q, V) = exp ( — f |V — AQ\ 2 1 is an invariant measure of (|2.26p . In contrast we do not know an 
analytical expression for stationary states of (|2.2ip for general choice of A. 

Second, we could consider the case of a non-homogeneous incompressible background flow, where the average 
velocity at the point q* in (|H1|) is u(q* ) where u is a divergence-free vector field. For any such u, the microscopic 
bath dynamics 

<iq = vdt, 

dv = Vu(q)vdi, 
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preserves any bath measure of the form (|2.16l) - (|2.17l) with v — Aq replaced by v — u(q). In this case, the 
dynamics of the large particle in the mechanical process converges, in the limit m — > 0, to 

dQ = Vdt, 
MdV = - 7 (V - u(Q))dt + adW. 

This extension likewise does not pose any theoretical difficulties, at least when u is Lipschitz continuous. This 
reduces to the case presented above when u(q) = Aq, where divu(q) = is equivalent to tr A = 0. 

Finally, a much more challenging extension is the case of multiple large particles. The restriction to a single 
large particle is a necessary assumption for the argument here and in [HEM?], since it allows us to estimate 
the distribution of velocities of the bath atoms that collide with the particle. In particular, fast moving bath 
atoms can collide at most once with the large particle, whereas in the multi-particle case these atoms could 
bounce between large particles and possibly collide with them many times. Extending the proof of Theorems 1 2. II 
and 12.21 to the case of multiple particles immersed in a bath is a non-trivial task, as atoms of any speed can 
bounce between the large particles, and there is also a shadowing effect that particles have on one another. 
Kotelenez |13j treated the multi-particle case using a mean field interaction representing a mesoscale interaction 
and showed how the bath can generate close-range forces among the large particles. Kusuoka and Liang treated 
the multi-particle case as well, using a weaker bath-particle interaction which prevents atoms from bouncing 
back and forth among the large particles [15] . We do not pursue a multi-particle derivation in this work; however, 
we do perform numerical tests of a natural extension of the limiting equations (|2.21[) to the multi-particle case 
in Section [3l 

3. Numerics for the multi-particle case 

We consider the application of the derived equations (|2.21j) to the simulation of a fluid composed of many 
identical large particles. We note that the derivation for Theorem 12.21 only applies to the case of a single 
particle, with no forcing except through the bath, whereas we now apply the equation to a system with many 
large particles that interact with one another. The many-particle case is the case of interest in applications, and 
we numerically investigate the qualitative behavior of the proposed dynamics to test its suitability for sampling 
flows with a given mean behavior. We consider a system in K 3 consisting of iV large particles whose position 
and velocity (Q, V) £ M. eN evolve according to the dynamics (|2.24[) which we recall here: 

dQ, = V. t dt, 
MdV, = -V Qi £(Q) dt - 7 (V S ; - AQi)dt + adW t . &' 1 ' 

The index i — 1, . . . , N runs over all particles, whereas the lack of index in the argument of E denotes the fact 
that it is a function of the full vector Q £ M. 3N . The Fokker-Planck equation for ()3.1I) is 



dtp = - Vq • (Vp) + M^Vv ■ (V£(Q)p) + Af^Vv • fr(V - AQ)p] + ^M~ 2 a 2 A v p. (3.2) 

We do not have an analytic expression for solutions of the time-dependent equation (|3.2[) or any steady-state 
solutions. Thus, it is useful to carry out numerical experiments on the response of the multi-particle system to 
the background forcing and explore the resulting constitutive relation. 

3.1. Simple Lennard- Jones fluid 

We now run numerical tests on a 3D flow of particles with Lennard- Jones interactions and background shear 
flow 

" s 

A= 



A closely related dynamics, differing in the choice of boundary conditions (we detail our choice below) and how 
the external forcing is handled, is considered in [12 where the authors perform rigorous asymptotic analysis on 
the invariant measure as well as numerical viscosity experiments. 
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We let fl = [-L/2, L/2] 3 be the computational domain, and we initially arrange N particles on a uniform 
cubic lattice with spacing a, where L = aN 1 ' 3 . The initial velocities are random, drawn from the measure 
Z^ 1 exp(— £|V — AQ| 2 ). We note that over the long times that we simulate, the solution is insensitive to 
the choice of initial conditions (we have also tested with zero initial velocities). The potential E denotes the 
Lennard- Jones potential energy with cutoff 

N N 

£(Q) = EE #»•«)> 



where 



4£ (^-^) +Cir + C2 ifr<r -*- 



Hr>r cut . 

The constants c\ and ci are chosen so that <fi(r) and 4>'{r) are continuous at r cut . The distance r,j = |Qi — Qj| 
is computed taking into account the boundary conditions which we now make precise. 

We apply the Lees-Edwards boundary conditions [l][9] to the system, which generalize standard periodic 
boundary conditions to be consistent with shear flow. At time t, the particle (Q(£), V(£)) has periodic replicas 

(QiW,g 2 (i),g 3 (i),V r 1 (t),T/ 2 (i),T/ 3 W) 

= (Qi(t) + mL + tnsL, Q 2 (t) + nL, Q 3 (t) + kL, V t (t) + nsL, V 2 {t),V 3 (t)) 

for all m,n,k € Z. This ensures that the periodic replicas of the system move consistently with the shear. 

We first perform a single run with the parameters e — 1,N — 1000, /3 = 1.0, s — 0.05, M = 1, o = 
0.7 _1 / 3 ,r cut = 2.6,7 = 1,0" = \/27/3 _1 . For this choice of parameters, the Lennard- Jones particles are in a 
fluid regime (see e.g. [22]). We run the simulation up to time T = 500, with a stepsize of At = 0.005 using the 
following splitting algorithm. We let (Q",V™) denote the approximate particle position and velocity at time 
t = nAt. We let a = e~* A *, and for i = 1, . . . ,N we define (Q" +1 , Vf +1 ) by 

v n+l/2 = y „ _ ^ M -l VQ ^ (Q n ) 
xn+1 _ n , a .-.rn+1/2 



q„+1 = Q« + AiV 

At 
~2 



V* = V^ +1/2 - ^Af- 1 V QiJ B(Q™+ 1 ) 



l-a^ 1/2 



Vr+ x - aV* + (1 - a )AQ™ +1 + (^^^ j GJ 

where G" ~ Af(0,I) where / denotes the dxd identity matrix. This is a composition of a Verlet step applied 
to the Hamiltonian portion 

dQi = Vidt, 
MdVi = -V Qt E(Q)dt 

with an exact integration of the remaining terms which represent the effect of the heat bath on the large 
particles, 

dQ l = 0, 
MdVi = -7(Vj - AQi)dt + adWi. 

3.2. Mean flow and stress 

In Figure [5] we display the results of our numerical simulation of shear flow. We partition the domain into 
K = 100 slices, 

_ L /2 + k —±L, -L/2 + |U 



R k = [-L/2,L/2]x 



K ' ' K 



x[-L/2,L/2] (3.3) 
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Figure 5. (a) The mean velocity (|3.4I) of the particles within horizontal slices (|3.3[) for 
a shear flow with strain rate s — 0.05, averaged in t. We plot 20 of the 100 bins used in 
creating the statistics, (b) The distance Q3.5P between the computed mean velocity and 
the background flow is plotted versus Q2 at the final time tg = T = 500. (c), (d), and (e) 
The variance of particle velocity around the mean is plotted versus Q2, where we again 
average in Qi 1 Qz 1 and t. This is well centered around /9 -1 . (f) The distance (|3.6|) of the 
computed mean velocity to the background flow is plotted versus time in a log-log axis. 
We observe the expected 0(£ -1 / 2 ) decay of the distance. 




10 100 1000 

Time t 



for k = 1, . . . , K. At any time ti = £At, we define a time average of the mean velocity in each slice k by 

v(u, k) = £"=oE^i v rWQD (34) 



eLoELwq?) 

where V" denotes the velocity of particle i at time nAt and where !#,. denotes the indicator function of the 
slice Rk ■ Since the flow is uniform in all but the e2-direction, averaging on slices helps improve the convergence 
of the flow statistics. The time average over the full simulation of the mean velocity is plotted in Figure [5ja), 
and we observe a linear profile equal to the applied background flow. 
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FIGURE 6. Shear stress ai 2 vs velocity gradient. We observe that the viscosity increases with 
the heat bath parameter 7, though for the two smaller values of 7, the viscosity is extremely 
close. For the values 7 = 0.1, 7 = 1.0, and 7 = 10, we find viscosities 77 = 1.2175, r\ = 1.2254, 
and 77 = 1.9518, respectively. 



We define the distance of the mean flow in each slice to the background flow, 



where Q 2 (k) = -L/2 



distfe k) = {VAU, k) - sQ 2 {k)f + V 2 2 (U, k) + Vffa, k) 



fc-1/2 



1/2 



(3.5) 



L is the y-coordinate of the center of rectangle R k - We plot the distance as a 



function of Q 2 in Figure [5jb) and we observe that it is uniform throughout the domain. The variance of the 
velocity is computed for each slice and is displayed in Figure EJc), (d), and (e). The variance closely matches 
the expected variance due to the background temperature, -r^-g = 1. In particular, it is statistically the same 
in each direction so that we observe a scalar temperature in contrast to the laminar flow case of (J2.13I) . In 
Figure Hff) , we plot the distance of the mean flow V(t^, k) to the background flow versus time, 



dist(if 



if: 

fc=i 



(y k ,i{t t ) - S Q k , 2 ) 2 + Vk,2 + 



1/2 



v. 



A-..-S 



(3.6) 



We observe an 0(i -1 / 2 ) decay in the distance, which is the convergence rate to the mean expected for Monte 
Carlo empirical averages. 

In Figure [51 we display the &i 2 term of the shear stress, computed using the Irving-Kirkwood formula |11) 
which we average over all time steps 



er = 



N 

E 



/ 



i N 

M(Vr - AQ?) (8 (V? - AQ1) + - "£(Qi ~ Qi) ® / 



\ 



«) 



3=1 



where 



/ 



(ij) 



IQi-QjD(Qi-Qj) 



iQi-Qil 

denotes the force on atom i due to atom j and |fi| = L 3 is the volume of 51. We note that it is standard to 
subtract the local average velocity from the pressure term |11) : hence we have subtracted the background flow 
at each step. We run ten independent simulations for each choice of the strain rate s = 0, 0.01, 0.02, . . . , 0.07 
and the parameter 7 = 0.1, 1.0, 10. The simulations are all run to time T = 500, with N, /3, M, a, r cut , and 
a chosen as above. We plot the stress along with error bars equal to 3 times the standard deviation. Time 
averages of the shear stress a\ 2 are plotted against the strain rate s. We notice a linear relation, and define 
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the shear viscosity r\ = —212._ Fitting the data, we find viscosity (with standard deviation) r\ = 1.2175 ± 0.0085 
for 7 = 0.1 whose confidence interval overlaps with the value r\ = 1.142 ± 0.087 reported in 22, Table IV]. 
The algorithm successfully simulates a system out of equilibrium, and the computed viscosity is consistent with 
previous computations for Lennard- Jones fluids. 

Appendix A. Convergence in the small mass limit 

In this section, we outline the proof of the convergence of the mechanical process (Q m , V m ) to the solution of 
the nonequilibrium Langevin dynamics (|2.21|) . as stated in Theorem l2.2l We construct a Markov approximation 
(Qm,V m ) that only counts the 'fast collisions' of the mechanical process (which we define below). The Markov 
approximation acts as an intermediate process between the mechanical process and the nonequilibrium Langevin 
dynamics. We prove in Appendix [Bl that the Markov approximation (Q m , V m ) converges to the solution of the 
nonequilibrium Langevin dynamics (|2. 21|) in the sense that 

(Qm, V m ){ <t<T} — >(Q, V){ <t<T} as m -* 0. (A.l) 

In Appendix [D] we show that the mechanical process and Markov approximation are close to each other, in the 
sense that 

(Qra,V,„)j <KT) - (Qm,V m ) {o <K T }A0 aS 771 -> 0. (A. 2) 

This is shown by coupling the Markov process to the mechanical process, that is for each realization of the 
mechanical process we associate a (random) set of realizations of the Markov process that are (with high 
probability) close to the realization of the mechanical process. By a standard result in probability theory [31 
Theorem 3.1], properties (IA.1I) and (|A.2|) allow us to conclude Theorem 12.21 that is, that 

(Qm,V m ) {0 < t < T} — KQ,V) {0 < t < T} as m ^ 0. 

The proof here is structured as in [7] , with changes made to handle the fact that the distribution of the initial 
bath atom velocities depend on space. We work with the pairing (Q m , V m ), since in our case, we cannot create 
a Markov approximation for the velocity alone. We give a self-contained proof, while also pointing out where it 
differs from the original argument in [TJ. 

A.l. Rate of fast collisions 

As noted in Remark 1 2. 11 the mechanical process is not a Markov process itself since a bath atom may collide 
more than once with the large particle or certain collisions may be impossible due to past motion of the particle. 
However, the expected value of the relative speed of a bath atom is OimT 1 ' 2 ), which is large since m is assumed 
small, whereas the expected value of the relative speed of the large particle is much smaller, 0(1). Most collisions 
happen between a fast-moving bath atom and a slower-moving large particle, and we show below that for such 
collisions there is no chance of recollision between the large particle and the bath atom involved. This motivates 
the following definition of "fast collisions" and the introduction below of a stopping time on the trajectory of 
the mechanical process that controls the large particle's velocity and position. 

Definition A.l. We call a collision a fast collision if the normal component of the relative velocity (|2. 19|) of 
the bath atom before the collision occurs satisfies 

v n > c m := m~ 1/5 

where v n = v • e n . Every other collision is called a slow collision. 

The particular scaling of c m is chosen so that c m — > oo and c^rn 1 ' 2 — > as m — > 0. We use the second limit 
in Appendix ID. 41 to bound the total effect of slow collisions. 

Definition A. 2. Let T > be given. For a given realization of the mechanical process (Q m , V m ), we define 
the stopping time 

r m =mmfinf{*:|L4|| |Q m (i)| + |V m (i)| > c m /8} ,Tj . (A.3) 
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FIGURE 7. The differential area carved out by fast atoms that collide with the large particle 
in the time interval [t,t + dt], see (IA.4[) . 

In the following, we assume that m is sufficiently small so that the initial condition satisfies ||A|| |Q(0)| + 
|V(0)| < c m /8, where we recall that the initial condition is independent of m. 

In Appendix [E] we show that almost surely r m > and that the mechanical process is well-posed up to r m . In 
particular we show that almost surely, before time r m there are only finitely many collisions and none of these 
collisions involve more than one bath atom at the same time. Intuitively, this means that the fast collisions 
have no memory of the large particle's trajectory, and furthermore it means that the rate of fast collisions is 
governed only by the initial bath configuration (J2.16I) . In Appendix [C] we show that if a bath atom experiences 
a fast collision with the mechanical process's large particle sometime in the interval [0,r m ), then this is the 
only collision that the bath atom undergoes in [0,r m ). From the coupling that we construct in Appendix [Pi 
between the mechanical process and the Markov approximation and a consideration of the expected size of the 
Markov approximation f (|D.9|) and (JD.10I0 . we can show that lim m _j.o P(Tm = T) = 1, so that for the sake of 
our convergence proof, on our time interval of interest, [0,T], fast collisions do not introduce memory terms. 
We thus define below a rate of collisions that the large particle experiences whenever the bath has configuration 
distributed according to (|2.16|) . 

Lemma A.l. Suppose that a large particle is at all times surrounded by a bath of atoms distributed according 
to (|2.16[) . Then the instantaneous probability of a collision between the particle at position Q and velocity V 
and a bath atom with velocity in the ball B(v; dv) on the surface element R d ~ 1 d^l centered around q = Q — Re n 
in the time interval [t, t + dt] is given by 

r m (v, e n , Q, V) dv dQ dt =A„ l i? rf ~ 1 max(u n - V n , 0)/ m (Q - Re n , v) dv dfl dt (A.4) 

where R is the radius of the large particle, \ m = m~ 1 < 2 A is the expected number of bath atoms per unit volume, 
v n = v • e n is the normal velocity of the incoming bath atom, and V n = V • e n is the normal velocity of the large 
particle. 

Remark A.l. The assumption that the large particle always sees the same environment is too strong, and indeed 
does not hold for the large particle in the mechanical process. However, as we show in Appendix [Cj for bath 
atoms that undergo fast collisions with the large particle, their pre-collision distribution does satisfy (|2.16|) . 
Thus, the rate (|A.4[) is correct for fast collisions of the mechanical process. 

Proof. In Figure [7| we sketch the differential volume element of size R d ~ 1 dH,(v n — V n )dt with base on the surface 
of the particle and height determined by the velocity of the incoming bath atom. The velocity of the bath atom 
and the large particle are constant on the time interval dt, so this surface gives the infinitesimal rate r m of fast 
collisions. □ 

A. 2. Rate of jumps and the Markov approximation 

We now use the rate (|A.4j) to define a rate for a Markov approximation to the mechanical process, where we 
only count the effect of fast collisions. We change variables to be in terms of jumps in velocity for the large 
particle rather than velocities of bath atoms. We first integrate out tangential directions v t from r m , and define 
the marginal 

f X {x) = / f(xe n +v t )dv t , 
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where we recall that / was assumed to be rotationally invariant, so that the left hand side does not depend on 
e n . We define also the scaled marginal, 

fa(x)=m 1 /*f x <jn 1 / 2 x). (A.5) 

Rewriting the collision rule (12.71) . we have that the change in normal velocity of the large particle in the en- 
direction, V = V^ — V n , due to a bath atom with normal velocity v n satisfies 

■~- 9m 

V = T7—( v n- V n)- ( A - 6 ) 

M + m 

We integrate out the tangential direction v t from the rate r m , make the change of variables from v n to V 
in (|A.6|) . and introduce a Heaviside function to restrict to jumps caused by fast collisions. The rate on jumps is 

f m (V, e n , Q, V) =A ro i? d - 1 (^P) 2 H {^^V + V n - (Aq) n - c m ) 

x max(F, 0)/i {-^V + V a - (Aq) n 

where 

q = Q - Re n 

represents the position on the large particle surface where the bath atom collides. Due to the Heaviside function, 
the minimum jump size V with nonzero probability is 

^ 2/77 

Knin = -jrjr- max ((c m - V a + (Aq) n ) , 0) , (A.8) 

M + m 
where we have kept implicit the dependence of V m i n on Q and V. We let 

A m (Q,V)=/ / r m (V,e n ,Q,V)dVdQ 
denote the jump intensity, and a short calculation shows that A m (Q,V) is independent of (Q, V) whenever 

\\A\\ |Q| + |V|<c m /8 

holds, due to the Heaviside function in f m . 

Definition A. 3. We define the Markov process (Q m ,V m ) starting from initial condition (Q m (0), V m (0)) = 
(Qm(0) ( V m (0)) to be a Pois son jump process where the velocity experiences jumps Ve n with intensity A m (0, 0) 
and with probability density junction 

A ro (0,0)- 1 f ro (F,e n ,Q m ,V m ), 

and between jumps, (Q,V) evolve according to (|2.6[) . 

Appendix B. Convergence of the Markov approximation to the SDE 

To show Theorem 12. 2\ we first show that the Markov approximation (Q m , V m ), with rate given by (|A.7|) . 
converges in law to (|2.21|) . We write the generator for the Markov approximation and show that it converges 
as m — > to the generator for (|2.21|) . These calculations are carried out explicitly below, allowing us to get 
the coefficients for the SDE in terms of the parameters of the heat bath and the large particle. The setup and 
calculations follow as in [7]. Because of the Q-dependence of the process, we work with test functions ?/>(Q, V) 
that are functions of the position and velocity. The inhomogeneity of the velocity field does not change the 
coefficients 7 and a compared to the zero background flow case (|2.9[) . 
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Lemma B.l. Let T > be given. The Markov process (Q m , V m ) satisfies 

(Qm,V m ){o<KT} KQ- V ){0<t<T} as m -> 

on [0, T], where (Q,V) satisfies the nonequilibrium Langevin dynamics (12.21)1 . 

For a stochastic process, we let T t : Co(R 2d ;R) — > Co(R 2d ;R.) denote the evolution semigroup, where 
Co(R 2d ;lR) denotes the set of continuous functions with zero limit at infinity. We define the infinitesimal 
generator L by 

^(Q,v) = to r «^Q,v)-mv) (B . 1} 

where the domain of L is the set of Co functions such that the above limit exists in Co- In the following, we use 
L m to denote the infinitesimal generator for the Markov approximation (Q m , V m ) and L for the nonequilibrium 
Langevin dynamics (|2.21|) . We prove Lemma IB. II with the use of Lemma IB. 21 below, which is a general lemma 
relating convergence of generators to convergence in law of the process. We recall that T>([0, T]) represents the 
space of cadlag functions on the interval [0, T]. A linear subspace K of the domain of L is called a core for L if 
L is the closure of L\k [13] . One can show that C^°, the set of compactly supported C°° functions, is a core 
for the nonequilibrium Langevin dynamics (|2.21[) using the ideas of [14) . 

The process (Q,V) is a Markov-Co-process, which means its evolution semigroup satisfies T t Co C Co and 
limt-K, ||IW> - VII = for all ib e C . 

Lemma B.2. Consider a family of Markov processes Z m on 2?([0,T]) with generators L m . Suppose that the 
Markov-Co -process Z has generator L. Let K be a core for L such that tp G K =>• rp G T>(L m ) for all sufficiently 
small m. Suppose that the initial distribution of Z m converges weakly to that of Z and that 

Vi/) G K, lim sup \L m ip(z) - Li/j(z)\ = 0. (B.2) 



Then 



Z m -^fZ m^O. D 



The above lemma is (7j Lemma 4.1], which can be found in similar form in |14) or |23) . 

B.l. Generator for the Markov process 

We apply the generator L m for the Markov process to a C^°(M. 2d ) test function, expand in powers of m, and 
show that we have, to leading order, a drift and diffusion term. 

Applying (IB.1[) to the Markov process, we have for any ip(Q, V) G C^°(R 2d ) that the generator satisfies 

L m ^(Q,V)=V-Vo>(Q,V)+ f f [^(Q,V + Ve n )-^(Q,V)}r m (V,e n ,Q,V)dVdn, 

where f m is defined in (|A.7[) and V m [ n , which we note depends on e n , is defined in (jA.8[) . We recall, from the 
scaling of bath atom velocities (I2.17|) . that E m (|v|) is 0(m -1 / 2 ). From ()A.6|) this scaling implies that E m (V') is 
0{m 1 ' 2 ). To find the limiting generator in the limit m — >• 0, we expand in powers of V and write 

7/HQ, V + Ve n ) - V(Q, V) =Ve n ■ V V ^(Q, V) + ~ V 2 {e n ® e n ) : V^V^Q, V) 

+ - V" 3 (e n ® e n <g) e n ) • 3 • V^V(Q, V + V*e n ), 
6 

for some V* G (0, V) and where -3- denotes the third-order contraction product, which we apply to the 3-tensors 
(e n ® e n ® e n ) and V^. We then write 

L m ^(Q, V) =V • V Q V(Q, V) + C m (h ■ V V ^(Q, V) + i h ■ V V ^(Q, V) + K m ) , (B.3) 
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where, recalling that q = Q Re n , we have 

h =m- 5 / 2 / r \> 2 e n f m f^±^t>+(V-^q) n ) ,/r,/<>. 



/2 = m~ 5/2 / i C ^ 3 (e n ® e n ) £ (^T^ + (V - Aq) n ) df>dft, 

m -5/2 f r°° ^ ~ fM + m~ \ - 

K m =—^ / ^ 4 (e n ®e n ®e n )-3-V^(Q,V + y*e n )/^ __!_y + (V - Aq) n dVdn, 



with coefficient 



>d-i (M + m 



C m = XR d - 1 [—5—) ■ ( B -4) 

Note that in our notation, we have suppressed the dependence of I\, I2, and TZ m on (Q, V). 



B.l.l. Remainder term 

We begin with the remainder term TZ m . Since the test function ip belongs to C£°, we can restrict to the 
case where Q and V are bounded, and we may assume m is sufficiently small so that |V| + ||^4|| |Q| < c m . 
This, along with the fmiteness of the moments (|2.3[) . allows us to estimate the order of lZ m . Letting x = 

(%$sV + toV2 ( V - M)n) , we compute 



\K m \<Cm- 2 

Cm 1 ! 2 f f x 



f r V 4 V^(Q,V + y*e„) f(M+!lv + m^ 2 (y-A q ) n )dVdn 

JS*-! JVmin V 2m ' J 

(7 rr ,l/2 f f°° -. 

<ttt—^e / (x-m^CV-AqW 4 V^V(Q,V + T/*e„) f 1 (x)dxdfl. 



We note that since we have assumed that m is sufficiently small, the minimum value of x is given by x m i r 
m x l 2 c m - Recalling that ip is compactly supported and that V* <E (0, V) = (0, ( jtt ) (x - m 1 / 2 (V - Aaf) n 



< 



we can bound Q , V in the innermost integrand to find the estimate {x— m ' (V— Aq) n ) VyVKQi V + T^*e n ) 

C(l + 2) where C may depend on ip but not on m. We find an upper bound on lZ m by extending the interval 
of integration and using the boundedncss of the marginals 



m l/2 l-oo 

l^-l ^ C 7^^¥ / f 1 + x ) f (*) da; ^ Cml/2 - ( B ' 5 ) 

(M + m) 5 J 



B.f.2. Diffusion coefficient 

We next turn to I2, and compute 



h = m- 5/2 f r (e n ® e n )V 3 /* f M±^F + (V - Aq) n ) d^df! 
= m- 2 I r (e„ ® ejF 3 / 1 (^4^t> + m'/ 2 (V - Aq) n ) dVdtl. 
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Let x — I 2^1/2 V + to 1 ' 2 (V — j4q) n j . We expand in powers of to, using the finiteness of the moments (|2.3 
and the fact that tp is compactly supported, giving 

i f\ r poo 

I 2 = 7TT-> u/ / (e n ®e n )(x-m 1/2 (V-AcL) n ) 3 f(x)dxdn 

(M + TOJ 4 J S d-iJ m i/2 Cm 
1 (a r /*oo 

= W LJ ( e n^^n)x 3 f 1 (x)dxdn + 0(m 1 / 2 ) 



= 174^ l (^ ® e ") d0 + ^(TO 1 / 2 ) 

16Sd_i,. r _, i/j, 
= ^d- $3l + 0(m/) ' 

where / denotes the d x d identity matrix and we recall that Sj-i denotes the surface area of the (d — 1)- 

1/2 

sphere S d_1 . The second line uses the estimate /J™ m x 3 f(x) dx < m 3 ^ 2 c 3 n J r /(x) dx = (^(m 1 / 2 ) along with 
a gathering of higher order terms in to. The third line uses the definition (|2.3[) of $3, and the fourth line uses 
the following identity 

/ (e n <gi e n ) dft = ~ I. 

Jgd-i d 

Multiplying by C m of (|B.4|) and letting to — > leads to the isotropic diffusion coefficient 

DJ = lim C m h = 4A ^" 1 f d - 1 $ 3 /- (B.6) 

In the case of /(w) = Z _1 exp ( — ^r ) with Z = 3^5, we have $3 = gs0-c and 

4V2XR d - l S d -! 



D 



^l 2 ^M 2 d 



B.1.3. Drift coefficient 

We now similarly expand I\ and find that the lowest order term, which is 0{m~ l < 2 ) 1 cancels out leaving an 
O(l) drift term. Indeed, we have 

h = j d i e„TO-t p V 2 f m (^^V + V n - (Aq) n ) dVdQ. 
As before, we let x = %£$V + m l / 2 {V - Aq) n . We obtain 

A= ° ,3 / e n / (x-m l / 2 (V-A q ) n ) 2 f l (x)dxdn 

f [Rm^ 1 / 2 r 00 1 fi r°° 

= / e « ^T7^ / x2 ft*) dx + Tk (-- T ( V - A ^ f 1 ^ dx 

J S d-i [ M A J m V2c m M d J 



dfl + 0{m 3/10 ), 



1/2 
where the error terms are dominated by f Q m xf(x)dx < m 1 ' 2 c m L /(x) dx = 0(m 3 ' 10 ). The first term 

vanishes, 

/ e„ / x 2 f 1 (x)dxdfl = 0, 

since J sd -i e n dfl — 0. We recall q = Q — i?e n , and let W = V AQ. We have the identities 



/ e n • We„ dQ = / (e n ® e n ) : W dtt = -^ W, 

/ e n (e n ■ e n ) dtf, = / e n dSl = 0. 

JS^-i JS^-i 
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Combining the above, we have 



/l = " 1 W (V " AQ)<i>1+O(m3/10) 



Multiplying by C m of (JB.4I) gives the drift coefficient 

-±(V-AQ)= lira C m h = - ^"'f "^ *i(V - AQ). (B.7) 

M ro->-0 Ma 

In the case of /(w) = Z _1 exp ( — ^- ] with Z = ^&, we have $1 = * and hence 

_ 2V2\R d - 1 S d - 1 

1 ' ~ y^d ' 

B.2. Stochastic limit 

Combining the expansion of L m in (|B.3|) . with the expressions for I\ and I2 (|B.6|) and (|B.7|) and with the 
bound on the remainder (|B.5[) . we have 



L m ^(Q, V) = V • VqV(Q, V) - M-^CV - AQ) ■ Vy^(Q, V) + \ DA v iJj(Q, V) + 0(m 3 / 10 ). 

We thus have a generator L m that in the limit m — > converges in the sense of (|B.2j) to the generator of the 
nonequilibrium Langevin dynamics (J2.21I) . where 7 is given in (|B.7[) and er = M\[~D is defined using (|B.6|) . in 
agreement with (J2.22J) . Thus, we use Lemma TB. 2 1 to conclude the convergence of the Markov process as stated 
in Lemma [B. II 



Appendix C. Fast collisions cannot lead to recollisions 

We show in this section that in the mechanical process until the stopping time (|A.3[) . any bath atom that 
undergoes a fast collision (recall Definition IA.1|) cannot recollide with the large particle and that the bath atom 
cannot have previously collided with the large particle. This shows the fast collisions experienced by the large 
particle have rate r m in (|A.4[) . This observation is necessary when coupling the Markov approximation to the 
mechanical process, which we do in Appendix [Dl 

Lemma C.l. For a given trajectory (Q m , V m ), suppose that the large particle experiences a fast collision with 
a bath atom at time t\ G [0,T m ). Then there are no other collisions between this atom and the large particle in 
the time interval [0,r m ), where r m denotes the stopping time (|A.3I) . 

Proof. We recall that the relative velocity of the bath atom is v = v — Aq. By the choice of bath dynamics (|2.20|) , 
the relative velocity only changes by colliding with the large particle whereas the velocity v changes with time 
according to (12.18[) . We consider a fast collision, 

v n > Cm = m~ 1/5 . 

We write the collision rule (|2.7[) in terms of relative velocity before making use of the bound on particle position 
and velocity in (|A.3|) (note that since the collision occurs on the particle's surface, q(ii) can be bounded in 
terms of Q(ii)). We then have the following bound on the relative normal velocity of the bath atom after the 
collision 

M-m 2M Tr 2M .. . ,, 

"n + T^ Vn-TT-, (-Aq(*l))n 



M +m M+m M+m 

2M tc 



< 


fM-m^ 

\M + mj 


U 


< 


2 

--Cm. 





M- 



R\\A\ 
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The last line holds for to sufficiently small. This shows that after colliding, the bath atom moves away from the 
particle with a large velocity. 

We let e n i denote the normal vector for the fast collision at time t\ and look at the future velocity of 
the bath atom, v(i) = Aq(t) + v'(i) for t £ [ti,T m ). From the above computation we see that the relative 
velocity of the bath atom is pointed away from the large particle in the e nl -direction. Due to the bound 
on the position of the large particle in (|A.3[) . any recollision with the bath atom must occur in the region 
11-^1 1 IQI ^ c m /8; however, whenever the bath atom is in this region, its velocity in the e ni i-direction satisfies 
v(i) • e n ,i < (c m /8 — 2c m /3) < — %»■. The large particle's velocity is bounded by c m /8, and is thus too low to 
overtake the bath atom. Therefore, there cannot be a recollision before time r m . Likewise, before the collision, 
the velocity in the e n> i direction satisfied v(i) • e nj i > 2&a ) which is faster than the large particle's speed which 
is bounded by c m /8, so it is impossible that there were previous collisions before time t = 0. □ 



Appendix D. Coupling the mechanical process to the Markov approximation 

In this section, we couple the Markov process (Q m ,V m ) with the rate (|A.7|) to the mechanical process 
(Qm,V m ), defined in Section |2~T1 That is, for each realization of the mechanical process, we associate a set of 
realizations of the Markov process that are, with high probability, close in the L°°-norm. For the coupling, we 
prove below that 

(Qm,V m ){ <KT) - (Qm, V TO ){ <t<T} >0 as TO ->• 0. 

As in |7] , we define a stopping time when the processes first differ by e > and bound the total effect of the 
velocity jumps not shared between the two processes to show that in the limit to — > the stopping time is 
greater than or equal to T with probability 1. 

As described below, for the majority of the velocity jumps in the mechanical process caused by fast collisions, 
we subject the Markov process to the same velocity jump. The construction of the Markov process differs slightly 
from that in [7]. Here, we couple jumps in the velocity rather than collisions with bath atoms, which makes the 
argument simpler. This simplification is made possible by the additional assumption that f(x) is decreasing, 
which we have assumed in our case in order to properly handle the fact that the bath atom velocity distribution 
depends on position. 

D.l. Coupling and convergence 

We construct the coupling up to the stopping time r m for the mechanical process (|A.3|) . and extend the 
definition of the Markov process up to time T, if necessary. Let / m0 ch = {£1,^2, ■•■} denote the set of all 
times up to r m at which the large particle in the mechanical process experiences a collision. This set is shown 
to be almost surely finite for any to in Appendix [El We define v n ^ = v n (ti), v n ^ = v n (ti), V n ,» = V n (ij), 
e n .j = e n (£,;), etc. We let J s iow = {ti £ /mcch : |v n ,i| < c m } denote the set of jump times due to slow collisions 
and /f as t = -Zmech\^siow denote the set of jump times due to fast collisions. These sets of jump times are random 
variables since the initial condition is random. 

For a given trajectory of the mechanical process (Q m ,V m ), we define (Q m (0), V m (0)) = (Q m (0), V m (0)), 
and for most of the fast collisions of the mechanical process we subject the Markov process to the same jump in 
velocity that the mechanical process undergoes. We selectively accept some subset of the fast collisions in the 
time interval [0, r m ) and add additional jumps in the velocity to ensure that the jumps of the Markov process 
have the rate f m (V, e n ,Q m ,V m ) defined in (|A.7[) . (Recall that the Markov process has no slow collisions.) 
More precisely, for every i; £ /fast, we choose to apply a jump with velocity change Vie nt i to the Markov process 
with probability 

PkBep(V^,e n) i,Q m ,V m ,Q m ,V m ) =min ^ ' -,1 . (D.l) 

For t £ [0, r m ), we add additional fast collisions to the Markov process with the Poisson rate 

r a dd(^,e n ,Q m ,V m ,Q m ,V m ) = maxff m (y,e„, Q m ,V m ) - f m (V,e n , Q m , V m ),0J . (D.2) 
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After accepting collisions with probability (|D.1[) . and adding new collisions with rate (ID.2|) . a short calculation 
shows that the rate function for the Markov process in the time interval [0, r m ) is f m (V, e n , Q m , V m ). If r m < T, 
we extend the Markov process to [0, T] by adding additional jumps with rate f m (V, e n , Q m , V m ). 

We let /Mark denote the set of all jump times of the Markov process, 7 rem C /fast denote the set of jump times 
of the mechanical process that were not chosen for the Markov process, and / a dd C /Mark the set of additional 
times at which jumps were added to the Markov process (in the whole interval [0,T]). From the construction 
above, we have that 

/Mark = (/fast \ /rem) U /add ■ 

We note that the sets of jump times /mech, /Mark, /slow, /rem, and / a dd are random variables. 

Since the chosen realization of the mechanical process may not be defined on the full time interval of interest 
[0,T], we make the convention that 



sup |V m (s) - V m (s)| = oo if T m < t. 
«e[o,t] 

In particular, part of the proof of convergence will be the fact that lim m _>.o P( T m = T) — 1. We then show the 
following convergence in probability result. 

Lemma D.l. For all T > 0, 

(Qm,V m )jo<KT) - (Qm,V m )j <KT) >0 as TO -¥ 0, 

where the convergence in probability is with respect to the L°°([0,T])-norm. Equivalently, for all T > and any 
£>0, 

lim P ( sup |Q m (i) - Q m (t)\ + \V m (t) - V m (t)| > e ) = 0. 

™^0 \ t6 [0,T] / 



To show Lemma ID. 11 we prove the following lemma which says that the interval of convergence can be 
extended by a finite time step. 

Lemma D.2. If to > is such that 



lim 1 

m— fO 



SUp Q m (t)-Qm(t) + V m (t)-V m (t) > 

vte[o,*o] 



for all e > 0, we then have 



lim I 

m->0 



for alle>0 where z = min [M(192\R d ~ 1 S d -i<i> 1 ) 



SUp Qm(<) ~ Qm(t) 
K te[0,t +z] 



V m (t)-V m (t) >e =0 



' 2(1 + ||A|| 



Since the initial conditions for the mechanical process and Markov process are the same, the hypothesis of 
Lemma [D.2I holds for to = 0. Thus, Lemma [D.l I follows immediately from Lemma [D. 21 and we are left with 
proving Lemma lD.21 The computation justifying the particular choice of z is performed in Appendix ID.3.41 

D.2. Error decomposition 

We prove Lemma ID. 21 by splitting the error in the velocity into two contributions and bounding them 
individually. We fix e > and to > 0, and in the following we denote by / mcch = /mech H [to, to + z], /^ark = 
/Mark H [to, to + z], and likewise for /g low , I^ cm , and I' add . For times t e [to,T m ), the error in velocity between the 
Markov and mechanical processes is 



V m (t)-V m (t) 



V m (t )+ J2 ^ e n,i l[to,t](*i)-V m (t ) 



^2 y^n,i l[to,i](*» 
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where we recall that U is the time that the jump Vie n! i occurs and where Is denotes the characteristic function 
of the set S. 

We split the error into two sources. The first source of error is the change in velocity due to slow collisions 
of the mechanical process, since the Markov process does not include any slow collisions. We denote this 
contribution in the interval [£q , t] by 



Wslow(t) 



Y^ ^ e *M l[to,t](*i) 



(D.3) 



The second source of error is the change in velocity due to the jumps that we added and removed when coupling 
the Markov process to the mechanical process. We dehne 



Wex(*)= J2 Vi t [t0 , t] (ti) + Y, ^ l[*o,t](*i)- 



(D.4) 



By definition, W cx (t) > for all t (recall that V > 0). The error terms Wslow and W CK are random variables. 
We decompose the error in the velocity as 



|V m (t)-V m (i)| = 



V m (t )+ Y V* e nA{t ,t]{U) - V m (*o) - Y ^ e n.» 1 [to,t](*«) 



tier 



t.eii 



< V m (t )-V m (t ) 



Y Vi e n,iMto,t](ti) 



U&li, 



(D.5) 



+ Y \Vi e nAlto,t](ti)\+ Y |^ e n,jl[to,*](**') 

|v m (to) - V m (t ) + W s i ow (t) + W ex (t). 



We fix e > and, as in [7 , define the stopping time 

T^minf inf if. V m (t)-V m (t) >e/2\,t + z) 



(D.6) 



The use of a stopping time gives us a bound on how close the mechanical and Markov processes are, which allows 



us to then bound the difference in jumps each process experiences. Provided that we have 
e/4, which we can assume from the hypothesis in Lemma TD. 21 we have that 



sup 

*6[*o,r* 



Qro(i) — Qm(i) < Qm(io) — Qm(*o) 



2 SUp 



Qm(*o) ~ Qrat'o, 

V m (t)-V m (i) <e/2, 



< 



where we recall z < 2(i+ \\ a \\ ) • Thus, we have the relations 



sup 
te[t ,T* 



Q m (t)-Qm(i) + V m (t)-V m (t) 



<e, 



sup || A|| Q m (t)-Q m {t) + V m {t)-V m (t) 
*e[* ,T^) 

From the above arguments, we may show Lemma ID. 21 by showing that 



< e. 



(D.7) 



lim ¥({ T ; n < t + z}) = 0. 

m— »0 



(D. 



2G 



It is difficult to estimate a priori the probability P(r m < to + z), which appears indirectly in the estimate of 
P(r* n < to + z) through the definition of T m (|D.6[) . To aid in the estimate we define the set of trajectories where 
the Markov process remains small compared to c m is 



G m ={ sup [\\A\\ Q m (t) + V m (t) <c m /16 

lO<t<T v ' 

From the convergence of the Markov process in Lemma IB. 11 we have that 

lim P(G m ) = 1. 

m— »0 

By the triangle inequality, we see that for trajectories belonging to G m , r m <to + z 
e/2. Thus, in order to have (|D.8j) . it is sufficient to show that 



(D.9) 

(D.10) 
V m (r*) - V m (r*) > 

(D.ll) 



(D.12) 



lim P({r m <t + z}n G m ) = 0. 

To show (|D.11[) , we use (|D.5|) to break apart the terms: 

{r* m <t + z}n G m c {{W slow (r^) > e/6} n G m } U {{W ex (r m ) > e/6} n G m } 
u{|v ro (to)-V m (t )| >e/6}. 
In terms of probability, we have 

P({< < *o + z} n G m ) < P({W s w(r m ) > e/6} n G m ) + P({W ex (T m ) > e/6} n G m ) 

+ p(|v m (to)-V m (f )|>e/6). 

By the hypothesis of Lemma ID. 21 

lim P ({ V m (t ) - V m (*o) > e/6J) = 0. 
We prove convergence of the other two terms in the following sections. 

D.3. Different collision rates 

Since we restrict our attention to times t < r m , Appendix [Cl shows that the rate of jumps due to fast collisions 
in the mechanical process is given by (IA.7|) . We now estimate 

W r «(T^)= J2 ^ 1 [*o,r*]fe)+ J2 ^IWlW. 

rem ' *- add 

which bounds the total effect of the jumps that were added and removed in the coupling process using the 
expressions in (|D.1[) and (|D.2|) . These jumps have the rate function 



1 * 7 ^11 7 ^lm 7 * 777 7 ^C77l 7 * 771 I ^TM I * 7 ®11 7 *** 



67777 * 771 



V m - rv 



(^,e„,Q 



V r . 
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The goal of this subsection is to show the following. 

Lemma D.3. Fix e > 0, and let r* n be given by (|D.6[) . Let G m be given by (|D.9[) . and let W ox be given by (|D.4[) . 
We have that 

lim P({W e x(r m ) > e/6} n G m ) = 0. 

771—^0 

Proof. In the following, we bound the error W oyi in terms of a pair of compound Poisson processes, whose 
definition we recall below. The error W ox is not a compound Poisson process itself, due to the dependence of r ex 
on Q m , V m , Qt77, V m , and e n , but we build a compound Poisson process from realizations of Wex, by increasing 
the size of certain jumps as well as adding additional jumps as detailed below. 
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D.3.1. Compound Poisson Processes 

We begin by recalling a few properties of compound Poisson processes before proving our bound on W ex . 
(Definitions may be found for instance in [5].) 

Definition D.l. A compound Poisson process J is a stochastic process defined in terms of its rate function 
r(V) by 

N(t) 
J(t) = J2 Vi 

i=\ 

where N(t) is a Poisson process with rate A :— f„r(V) dV that is independent from the jumps {Vi]^L l which 
are independent, identically distributed random variables with probability density function k~ l r(V). 

The expected value of a compound Poisson process satisfies 

E( J(t)) = E(N(t))E(V) = t J Vr(V) dV dt. 

Jr 

Using the independence of {Vi}i<i<N a ^ and N CXl we can simplify the variance 

'iV(t) 



E((J(t) - E( J(t))) 2 ) =E\[j2Vi- E(J(t)) j 



E | [ J2 ft - E(N(t))E(V) 

f N(t) \ 2 \ 

= E | [ J2 (Vi ~ HV)) + [N(t) - E(N(t))]E(V) 

/ =1 / / 

E \^2 & ~ E(U)) 2 + E ((N(t) - E(N(t)) 2 ) E(Vf 

E(N(t))E([V - E<y)] 2 ) + E(N(t))E(V) 2 
E(N(t))E(V 2 ) 

t f V 2 r{V)dV. 



(D.13) 



To go from the third to the fourth line, we have used independence of N(t) and Vi, and to go from the fourth 
to the fifth we have used the independence of the jumps {Vi} i=1 and the following property of the variance of 
exponential distributions: 

E ((N(t) - E(N(t)) 2 ) = E(N(t)). 

D.3.2. Building the auxiliary process 

Since there are Heaviside functions within the definition of r mi we consider three regions: one where both 
rate functions are nonzero, a second where one rate function is identically zero, and a third where they are both 
identically zero. For t € [to,T^), we define 

vi = --j— — min (\c m -V n + (Aq) n ) , (c m -V n + (Aq) n )j 
v 2 = -rr— — max (\c m -V n + (Aq) n j , (c m -V a + (Aq) n )j 
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where we recall that q = Q m — Re n and q = Q m — Re n are the locations of the collision on the surface of 



"'" " " < «i < V 2 < -MT^t from <ED and CO), and v 2 - Vl < 



the respective atoms. Note that M+m — 
from (|D.7p . 

In the range t>i < V^ < W2 , we can directly bound r ox using the monotonicity of / 



Ira 

M+m c 



<(^,e n ,Q m ,V m ,Q m ,V m ) = C m m- 5 / 2 y/i(^-^(F- Ul ) + c ?) 

2m 

M + m 



< C m m- 5 / 2 Vfl ( 



V - — 
2m ' 8 



where 

d-i/ M + m 

We integrate out e n , giving the factor Sj-i, and define the rate 



L-m — Art 



,-, /• c -5/2x7,1 (' M + m g. c„ 

r x (V) = C m S d -im 0/ V/ m ( — V" - — 



To construct the Poisson process, we add additional jumps in V with the (positive) rate 



r i(V) - / r ex (F,e n , Q m ,V m ,Q m ,V m ) dfi. 



For v 2 < y < cxd, we have the following estimate: 



= a„m- 5 / 2 y 



/i ( ^? + K - (Aq) n ) -4 f ^V + V n - (Aq) n 



2m 



V 2to 



= C m m- 5 / 2 y 
< C m m- 5 / 2 V 



t i I M + m .- \ x (M + m ~ 

/m ( o___ (X - W 2J + C m -/, m — (V - Vl) + C„ 



2 m 



2to 



,! .M + rn.^ , \ ,i (M + m.~ . 



2m 



2m 



where the inequality uses the monotonicity of / . We cannot directly eliminate V2 from the argument by relying 
on monotonicity arguments since we have a difference of /, so we shift by replacing jumps generated by the 
process corresponding to r cx by larger jumps. For any jump of size v 2 < V < oo, we replace it with a new, 
larger jump of size V = V + M + n m g L ~ v i- We also increase the rate of jumps by a factor of X and integrate 
out e n , giving the rate 

r 2 (V) = CmSd^m-^V 



i M + m - Cm 
f 1 — — v — 

■'"' ! 2m 8 



{ M + m c^ 

Jm[ 2m 8 



The limits of the regions above, v\ and v%, also depend on (Q m ,V m ) and (Q m ,V m ). To remove this 



dependence we extend the interval of definition for the two rate functions defined above to 



M+m 8 — V — 

u — — iui , i vim uv^ M — - < V < oo for r 2 . This can be done by adding additional jumps to the auxiliary 
Poisson processes above. All told, we have build realizations of two compound Poisson processes 



2m 9c„ 



for ri and to 



Ira 7c„ 



Jm,k — 



iV exfc 

E 

i=l 



Vi l[t„,to+-](*0 f or k = 1,2, 
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with the respective rate functions 



f cx AV) = n(V) 1^,2^] (^±^V 1 



2?7J 



f ex , 2 (y) = r 2 (f)l [Z ^ r0o) (^^ 



where we have used the dummy variable V for both processes and where we suppress dependence on time in 
the notation J m .k- 

For trajectories in G m , the Poisson processes J m .i and J TOj 2 provide the bound 

Wex(Tm) < Jm,l + Jm,2- 

Therefore, 

P({W«(t£) > e/6} n G m ) < P({J m ,i > e/12} n G m ) + P({J m , 2 > e/12} n G m ). 

To show that P({W ex « l ) > e/6} n G m ) -> in (|DTT2]> . we show below that P({ J m , fe > e/12} n G m ) -> for 
fc = l,2. 

D.3.3. Bounding J m ,i 

Applying Markov's inequality, we have 

P({ J m ,i > e/12} n G m ) < yE m (4,i). 

We make the substitution x — ^ I+ JJk V — m 1 / 2 ^ and bound the integrand on the finite interval to obtain the 
bound 

/>00 

EmGW) = ((*0 + ^) - «0) / W ox .i(y) ^ 

Jo 

r^ c -5/2 r^^Mfi ( M + m - r Cm\ - 

= C m S d . x zm I J ^ 7 ^ Vf m [—-V- — \dV 

M+m 8 V ' 



< 



8G m 5 d _ 1 zm 1 /2 p ^ / x ^ 
§ml/2cm UV2 g J 



(M+m)' 1 ji m i/2 Cm \m L ' z 8 / 

8G m 5 , d _izm 1 / 2 ^9 Cl7 
(M + mf 

XR ~ Sd- 
32(M + m) 



81XR Sd-iz 1/2 2 



where we recall the scaling of f m (IA.5I) . We recall that c m = m 1 / 5 so that for any finite z we have 



Thus, we conclude 



D.3.4. Bounding J m ,2 
We now wish to show 



lim E m (J m ,i) = 0. 

m— fO 



lim P({ J m ,i > e/12} H G m ) = 0. 

m— >0 



lim P({ J m , 2 > e/12} n G m ) = 0. 

m— J-0 



:-'»() 



For this estimate, we first show the bound ^ m (Jm,2) < e/24 before invoking Chebyshev's inequality. We have 
that 



E m (J m , 2 ) = */ Vf„j(V)dV 

Jo 



C m Sd-izm 
8C m Sd-i 

2; 

(M + m) 3 



-5/2 / yi 

2m 7a 



,V2 



M + m 8 



II (^V - ^) - fl (^V - %r - -: 



V 2m 



2m 



1m 1 / 



i/a„ Vm 1 /^ 8 / 



fir 



^(fc^+e) 



(^ + T" £ ) 2/(a;)dX 



8C m Sd- 
(M + m) 3 



1 1/2 

— zm ' 



2£ 



^72 +f )-0/(*)*» 



Im. 1 /^ 



(s^ + t)'^w* 



m 1/2 (|c m +e) 

Vm 



4Afl rf ~ 1 S d _ 1 
M 



Z£$i + 0(c^m 1/2 ) 



dV 



(D.14) 



(^ + f" £ ) 2 ^ 3 



fix 



S^y - ^ and -?™ = M±^V - So. + e . We recall 

2m 8 m 1 / 2 2m S 



where we made two changes of variables, w, 

linim^o c m m 2 — 0, so that only the leading order term remains in the limit. 

Now, by our choice of z = min ( M(192Ai? d ~ 1 S' ( j-i < ! > i)~ 1 , 2(1 -nrjm ) in Lemma ID. 21 we have 



^ m (J m ,2) < e/24 



(D.15) 



whenever m is sufficiently small. We later need a lower bound on the expected value, E m (J mj 2) > Ce, where it 
is not necessary to precisely find the constant so long as C does not depend on e or m. Such a bound follows 
from the above equations and the definition of z. 

Applying Chebyshev's inequality and using (|D.f5[) . we have 



P({ J m ,2 > e/12} n G m ) < P({J m , 2 > e/12}) 

< P({Jm,2 > 2E m (J m , 2 )}) 
= P({J m ,2 - E m (J m , 2 ) > E m (J m , 2 )}) 
E m ((J m> 2 — E m (J m> 2)) ) 



< 



Em(Jm, 2 ) 2 
Em(A^ex,2)E m (y 2 ) 
Em(Jm, 2 ) 2 : 



where the final inequality follows from (|D.13|) . 
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We can then estimate E m (7V eX] 2)E m (V 2 ) with a similar computation as for E m (J TOi 2) 



E m (7V cx , 2 )E m (y 2 ) 

'to 

16C m Sd-i 
(M + m) 4 ' 



to-\-z poo 



V 2 f c ^ 2 {V)dVdt 



.3/2 






l/2f3 



(lc m +e) 



V TO l/2 



+ — - e f{x)dx 



(D.16) 



< 



12-S d _i 



M 2 



izm 3 / 2 / e (_|_ + ^L) /(^dx + oCm 1 / 2 ) 



< Czm 1/2 + o{m 1/2 ). 

We see that this converges to zero as m — >• 0. 

Combining the above, for m sufficiently small, we have the following estimate on J m ,2, 



P({ J m ,2 > e/12} n G m ) < 
and so in view of (|D.14[) and (jD.16|) . 



t (JVcx, 2 )E m (T/ 2 ) 
E m (J ro , 2 ) 2 



lim P({ J m , 2 > e/12} n G m ) = 0. 

Combining the estimates of this section, we conclude 

limP({We X «J > e/6} n G m ) < lim P({ J„ hl > e/12} n G m ) + lim P({J m , 2 > e/12} n G m ) = 0, 

m— >0 m— fO m— fO 

which completes the proof of Lemma ID.3I 



□ 



D.4. Slow collisions 



In this section, we bound the effect of jumps of the mechanical process caused by slow collisions. To do so, we 
track the total change in the relative velocity of a bath atom, v = v — Aq. A bath atom that undergoes a slow 
collision may collide with the large particle many times; however, if the normal component of the bath atom's 
velocity ever exceeds c m after a collision, then no further collisions are possible, which follows from arguments 
used in Appendix [Cj 

Suppose a particular bath atom undergoes k collisions with the large particle. We denote by v, (respectively 
v£) for i = 1, . . . , k the relative velocity of the bath atom before (respectively after) the ith collision. We show 
below that the total change in velocity for a bath atom is bounded by \v' k — Vi| < 8c m . By the preservation of 
linear momentum for an elastic collision, the change in velocity for the large particle due to all the collisions 
with this particular bath atom is V — V = S(vj^. — Vi). We also show below that the number of unique bath 
atoms experiencing a slow collision with the large particle in a given finite time interval is bounded above by 
Cm~ 1 / 2 c m - Thus, the total expected change in the large particle's velocity is bounded by Cm 1 / 2 ^ which 
converges to zero as m — >• 0. 

Lemma D.4. For r* n be given by (|D.6j) . G m given by (|D.9j) . and W s \avi defined above in (|D.3|) . we have that 

lim P({VK slow «) > e/6} n G m ) = 0. 

m— >0 

Proof. We focus now on a single bath atom and consider the set of times I roc = {i;}i=i....,fc of all collisions of 
this single atom in the interval [to; r m] (t ne f ac t that there are almost surely only finitely many collisions is 
shown in Appendix [E|. We later estimate the total number of distinct bath atoms that undergo slow collisions 
with the large particle in the time interval. We denote by e ni i the normal direction for the collision at time 
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Figure 8. A slow collision with high tangential velocity in the plane defined by the bath 
atom's tangential and normal velocities, (a) If the tangential velocity of the first collision is 
large all subsequent collisions must occur in one hemisphere, in the positive e^i-direction. (b) 
The later collisions have normal velocity smaller and smaller, since the large component in the 
e t .i direction is outward pointing from the surface of the particle at the site of subsequent 
collisions. 

ti. We denote by v n ,i: v' ni the normal velocity for the particle before and after the collision at time ti, and 
similarly for the tangential velocity. Finally, we let v n ^ = v n) j • e n j. 

We recall that in a slow collision, the normal velocity of the bath atom satisfies v n < c m , whereas the 
tangential velocity cannot be bounded, so we do not have a bound on |v|. In the following, we consider three 
cases, based on the size of the tangential velocity of the bath atom. 

For the first case, we consider that |v t] i| < c m for alH = 1, . . . , k. The sum of the change in v telescopes 
since v' n i - v ai = v^ - v; = v i+ i - v i; giving 



5X< 



£*«-*< 



= |vl-v 1 |<2V2c ?) 



As a second case, suppose that |v t ,i| > c m . Even with a large tangential velocity, it is possible that there 
is a recollision with the large particle due to the curvature of the particle (see Figure [8]) . We split the large 
particle into halves, with the equator perpendicular to et t as in Figure El The bath atom passes the equator 
and due to its high tangential velocity, which is unchanged by the collision, all subsequent collisions must occur 
on the same hemisphere (the shaded hemisphere in Figure Ufa)). We recall that the speed of the large particle 
is bounded by c m /8 (|A.3|) . Between subsequent collisions, the bath atom moves further in the e tj i-direction, so 
that (e ti i, e n .j+i) < (et,i,e n) j) and therefore (e t .i,e, v ;) < for all i. This in turn implies monotonicity of the 
normal velocity for the bath atom's collisions, 



V n ,i+1 



<v' 



(D.17) 



From Lemma lC-H any bath atom involved in a fast collision can never undergo slow collisions, so in particular 
we have that |v n i| < c m for i — 1, . . . , k. We use the monotonicity of the normal velocities (|D.17j) to form a 
telescopic sum on the normal component of the change in velocity, which becomes 



k 



A— 1 



< 


£ |Vn,t+l 
i— 1 


- V n! 


i\ + Wn : k' 


< 


V n ,l - «n,k 


+ K 


,k~ Vn,fe| 


< 


2c m + 2(c m + Cm/16) < 5c 



V,i,fe 
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To prove the third inequality, we use the bound on \v ni i\ < c m for i — 0, . . . , k, along with the bound on the 
large particle in (|A.3|) . and apply the collision rule (J2.7I) . 

For the third and final case, consider the case that after a certain number of collisions, |vt,i| > c m whereas 
until that time, we have that |vt,j| < c m and |v n ,i| < c m , for i — !,...,£— 1. This combines the previous two 



:W 



cases. The change in velocity satisfies 

k 



E V >M ~ Vn -' 



|V^_i - Vi| 



E v n,* - Vn,, 



E^ 



< 2V2c m + 5c m < 8c m . 

We estimate |vi| < V2c m and |v^_ 1 | < \[2c m as both normal and tangential components of the relative velocity 
are bounded. We bound the final term, involving collisions after the fast tangential velocity using the previous 
argument (|D.18j) . 

Taking the maximum over the above three estimates, we conclude that the total momentum change for the 
large particle due to collisions with this single, slow particle satisfies 



E M ( V .'m 


- v n)i ) 


= 


k 

E 


i=\ 






i=\ 



K, 



< 8mc„ 
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We now define N s \ ow as the number of first-time collisions that are slow collisions, and we employ a double- 
index notation for the collisions, so that the jth atom's ith collision has relative normal velocity v nj jj where 
1 < j < -Wsiow an( i 1 — * — kj- While the slow collisions are not a Poisson jump process, we note (|A.4| acts as an 
upper bound for the first-time slow collisions. (We recall that in the mechanical process certain slow collisions 
are made impossible by the past history, though no first-time slow collisions are ignored in (|A.4[1 .) Then we 
may estimate the number of distinct bath atoms involved in slow collisions by 



A 



E m (7V slow ) < z I I —-Bt-Hvn - K)+/„((v - ^q)„) dv n dfl 

J S d-i J m l ' z 

< XR^S^m-^z f m 2c m /i((v - Aq) n ) dv n 



< 2\R d ~ l S d _ 1 mr l l 2 



zc 71 



Combining (|D.19[) and (ID.20[) . we finally have the estimate 



P(Wiow(C)>e/6}nG m )< 




N slovr ki 

EE'i 

j=l i=l 



i,j v n,i,j 



>e/6 



P (^^N slow > e/6 



< 



< 



V M 
48toc. 



E m (A^ s i ow ) 



96Ai? d - i S' d _izc^m 



2 m i/2 



as m — > 0. 



(D.20) 



□ 



D.5. Conclusion of the proof 



Combining the above estimates on W s iow and W cx , we conclude that the left hand side of (|D.12j) converges to 
zero in the limit as m — > 0. This shows that Lemma lD.21 holds, and we then conclude that Lemma ID. II follows. 
Over finite time intervals, we can then conclude that (Q m ,Vm){o<KT} ~ (Qm,V m )jo<j<T) — ^0 as m — >• 0. 
Combined with the convergence result Lemma IB.1( we can therefore conclude the proof of our main result, 
Theorem l2~2l 
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Appendix E. The mechanical process is well-posed 

In this section, wc show that the mechanical process is almost surely well-posed for all times t < r m , where 
r m is the stopping time defined in (|A.3|) . Particularly, we show that as long as the large particle stays bounded, 
the degenerate cases of having a multicollision or having infinitely many collisions in finite time both occur with 
zero probability. Wc further show that almost surely r m > 0. We note that this well-posedness result is slightly 
weaker than the corresponding result in [6][7] for a heat bath with zero background flow since that result does 
not require a stopping time and is valid over finite time intervals that are uniform with respect to the initial 
condition. In contrast, the stopping time r m depends on the initial condition. 

We first show that whenever the particle is restricted to a ball |Q m | < £ for £ € N, then there are almost 
surely only finitely many atoms that can possibly interact with the particle in a given time interval [0, T], since 
only finitely many enter the ball. It is therefore sufficient to study the well-posedness of the dynamics for finite 
systems. After restricting to a finite system, we show that the mechanical process is almost surely well defined. 
Since |Q m | < c m /8 for t £ [0, r m ), we conclude that the mechanical process is well-posed at least until the 
stopping time r m . 

E.l. Finitely many atoms within a finite radius 

In the following, we work with fixed mass m for the bath atoms, which we absorb into other constants. 

Lemma E.l. For any finite time interval [0,T] and radius £ € N, there are almost surely only finitely many 
bath atoms that enter the ball |q| < £ in the time interval [0,T]. 

Proof. For a given initial configuration, let N k be the number of bath atoms such that k < |q| < k + 1 and 
|v| > sk, where we choose e — ±e~"' A '' +i z' T . We then have 

Em (N k )= / A m f m (q, v + Aq) dv dq 

Jk<\q\<k+1 J\v\>ek 

f lvl d+1 

= C[(k + l) d k d ]S d ^ / L\f{ m ^v) dv 

J\v\>ek l v l 

^rvmi lv| rf+ 7(v)dv 

(ek) d+i J M > rn i/2 ek 



< 



C_ 
k 2 ' 



where we note that C depends on e, but not on k. Then, almost surely X^fcLi ^k < °°- Suppose an atom 
has initial coordinates (qo,v ) where k < |qo| < k + 1 and |v | < ek for some k £ N. We have from the 
dynamics (I2.20P that 



d fh\ 2 \ . / , 



dtK"r)=« TM+ ^ 

|q| 2 -^|q| 2 -^|v 



> - 



2 1 I |2 1 1st |2 



This implies 



I^I^O^ + ^n^V-^'^-WS 1 ' 



2||aii+i; 2pn + i' 

>| qo |2 e -(2||A|| + l) t _| Vo |2 
> fc 2 e -(2||A||+l)t_ £ 2 fc 2 

= 3 fc 2 e -(2||A||+l)T 

4 
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Let £ £ N be fixed. We choose k^ sufficiently large such that -^fcoe^ ( -" j4 " + 2) T > I. Then, the only atoms that 
can enter the ball |q| < I must either be one of the finitely many such that |qo| < fco or be such that there 
exists k > ko such that k < |qo| < k + 1 and |vo| > ek. □ 

In a finite mechanical system, the first collision occurs at a time to > 0, since at time t — 0, no bath 
atoms are touching the large particle. Upon choosing £ = ^f-, this implies that the stopping time r m (|A.3|) for 
the mechanical process satisfies r m > 0. (Recall that that after Definition (IA.3I) . we have assumed that m is 
sufficiently small so that the initial conditions lie within the box that defines the stopping criteria. 

E.2. A finite system is well-posed 

In the following, we consider a system composed of the large particle along with N £ N bath atoms, in the 
finite time interval [0, r m ). We show that there arc almost surely only finitely many collisions, none of which are 
multicollisions. For any initial condition of the mechanical process we define the process on the time interval 
[0, r*) where r* is the first occurrence of either a multicollision or an accumulation of infinitely many collisions, 
or the first time the large particle leaves the box used in the definition of r m (see (|A.3[) V By definition, the 
mechanical process is well-posed on the time interval [0,r*), and we have t* < r m . From Lemma TE. 11 we can 
conclude that t* > 0, and we show below that almost surely t* = r m . 

We note that because of the bath atom dynamics (|2.18p , the kinetic energy of the mechanical process is not 
preserved. When restricting to a finite system, we can control the growth of the energy. 

Lemma E.2. In the time interval [0,t*), the kinetic energy of the mechanical system at time t, 

1 1 N 

K{t) = -M\(t) T V(t) + ^ mVj (t) T Vj (i), 



j=i 



satisfies K(t) < K(0) 



2\\A\\t 



Proof During an elastic collision, the kinetic energy is unchanged. Between collisions, we have 

±K(t) = MV(tf ±V(t) + £>v J -(t) T i Vi (t) 



N 

= J2mv j (t) T Av j (t), 

3=1 

and we thus have \-^K(t)\ < 2\\A\\K(t). Integrating gives the desired result. □ 

As a corollary, we have that for every initial condition, the speed v max = I — — ) bounds the 

maximum speed of any bath atom: 

max max |v»(*)j < v max . (E.l) 

te[0,T*)j'=i,...,JV 

We let Bn{£,K) denote the set of initial conditions where r* < r m , that is, those that lead to either a 
multicollision or infinitely many collisions, where N is the number of bath atoms in the system, where the large 
particle is restricted to |Q| < £, and where the total kinetic energy satisfies supr T i K(t) < K. If we show that 
the set Bn{£,K) has zero measure, then the set of initial conditions leading to such 'bad' collisions, which is 
given by 

OO CO oo 

|J |J |J B N (£,K), 

n=i e=i k=i 

has zero measure. 

We first restrict ourselves to the case of a system with the large particle and a single bath atom. 



3G 



Lemma E.3. For all initial conditions in Bi(£,K), there are infinitely many collisions in [0, r*) and the normal 
component of the atom and particle velocities before the ith collision satisfies 

lim \v ni - V b a\ = 0. 

l— >oo 

Proof. Since there is only one bath atom, the only possible bad collision is the accumulation of infinitely many 
collisions. Suppose a collision occurs with relative velocity v nj i — V n j = e, for some e > 0. Immediately after 
the collision v' n i — V£ i = —e. Since the system is composed of only a single particle and a single bath atom, in 
order for another collision to occur, the bath atom must accelerate until (v(t) — V(t)) • e n i > 0. 
From the bound on the bath atom velocity in (jE.ll) , 



— v(t)=Av< ||A||u max . 

We conclude that there exists S > such that the time between collisions satisfies ti+i — U > Se. Thus, an 
accumulation of collisions implies that lim^oo \v a ,i — Ki,i| =0. D 

Lemma E.4. The set Bi(£,K) has measure zero. 

Proof. For any e > 0, let t* denote the time of the first collision where \v n .i — V n: i\ < e. For P E N, let 
T^ C -Bi(£, K) be the set of initial conditions where t* € [pT m , ^^-T m ), for p = 0, . . . , P — 1. By Lemma TE.3[ 
an accumulation of collisions must lead to slow collisions, so that Bi(£,K) C UpJ^r^p. 

Up to time t* , the dynamics is well-defined, and the flow map (fit for the mechanical process preserves 
Lebesgue measure. For p = 1, . . . ,P - 1, 4>iL Tm (r£ p ) C r^ , so that |r£ p | < |r^ |. 

For initial conditions in r^ , we may suppose that P^ 1 is sufficiently small so that r* is the time of the 
first collision (recall from Lemma lE.31 that a collision with \v n> i — V n ^\ > e leads to a minimum time between 
collisions ti+i —W> 5e). Mapping back from a collision at time r*, we see that the bath atom's initial position 
satisfies |q — Q| < " m p T "' and the initial velocity satisfies \v n ^ — V n _i\ < ee" j4 " _ p L t; max . Because of the bounds 

on Q and K, we then conclude that \Tf \ < ^^e~p , where Cx,C% € R are independent of e and P. Summing 

over p = 0, . . . , P — 1, we have 

p-i 



\B x {t,K)\<\m l Km 

e— >0 P—yoo 



Uri 



p 

e,p 
p=0 



0. D 



Lemma E.5. The subset M C Bn(£,K) of initial conditions where the first bad collision is a multicollision 
has measure zero. 

Proof. Intuitively, the set of all multicollisions at a given instant r* has phase space co-dimension greater than 
or equal to two, so that the set of initial conditions leading to a multicollision has zero measure. We note that 
a multicollision could involve one or more bath atoms that have an accumulation of infinitely many collisions 
at time r*. 

For P € N, we let r^_ be the set of initial conditions such that t* G [-pT m , Pjr-T m ), where p = 0, . . . , P — 1, 
where we recall that t* is the time of the first multicollision. We note that since the flow map of the mechanical 
process preserves Lebesgue measure, |r^ | < |r^ |, so we only need to bound the measure of r^ . 

For initial conditions in r^ , two or more bath atoms must satisfy |q — Q| < t ' ma ^ T "' , so that we have 

|r£ | < w- We thus have 

p-i 



\M\ = lim 

P— foe 



p 

e,p 

p=0 



= 0. D 



Lemma E.6. The set of initial conditions leading to a bad collision has measure zero. 

Proof. We first consider the full set Bm{£, K). If the first bad collision in such a system is not a multicollision, 
then there is a time 5 such that in the interval [r* — 8, r*), the only collisions that the large particle experiences 
are with the bath atom that experiences the bad collision at r* . In this case, the arguments from Lemma IE. 4 
show that the set of such initial conditions has measure zero. Therefore, using Lemma lE.51 we conclude that 
Bn(£, K) has measure zero as well. 
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Wc then conclude that 



and that almost surely r* 



oo oo oo 



UUU b n{^k) 



= 0, 



r m and so the mechanical process is well-posed until at least the stopping time 

□ 



Appendix F. Convergence of the generator for the laminar heat bath 

In this section, we compute the limiting SDE for the two laminar heat bath models defined in Section 12.31 
The computations parallel those of Appendix [B] and we omit some of the details found there. The first heat 
bath model has velocity distribution of the form (J2.11I) , and the second has distribution of the form (|2.14p where 
we consider the specific shear flow 

" s 

A= 



for some s£K. We do not carry out a full convergence proof as we do for the model proposed in Section [ 
although such a proof could likely be carried out based on the same arguments as above — since the limiting 
dynamics obtained by the following formal computations has undesirable traits for sampling nonequilibrium 
states. 



F.l. Unidirectional flow 



We start with a background bath whose atoms are distributed initially according to (|2.1ip . where all bath 
atom velocities at initial time are oriented in the ei-direction. In fact, rather than working with (|2.11[) and its 
Gaussian distribution, we work with 



dfj, m (<l,vi) = X m f m (vi - sq 2 )dqdv 1 , 



(F.l) 



where we have the scaling f m (v) = m 1 ' 2 f(m}' 2 v), where the probability density function f{v) is assumed to 
be decreasing in v with four finite moments (|2.3p , and where we recall q 2 = Q 2 — -Re„ • e 2 ■ Note that we have 
eliminated v 2 and 1*3 in the above expression since they are identically zero. We recover (|2.11[) upon choosing 

f(x) = Z- 1 exp(-§x 2 ). 

The normal velocity for a bath atom before a collision is v n — v\ei ■ e n , and we keep the distinction between 
slow and fast collisions of Definition IA.1I The rate of fast collisions experienced by the particle is 

r m (vi,e n , Q, V) dvi dtt dt =\ m R 2 H((v x ei - sq 2 ) ■ e n - c m ) 

max(w n — V ni 0)/ m (i>i — sq 2 ) dv\ dQ, dt. 



In contrast to (|A.4[) . the rate function here includes a Heaviside function to ignore slow collisions. We let 
V = V„ — V n , which satisfies, in view of (|2.7|) 



V 



1m 



V m 



M + m 

1m 

M + m 



(vi ei -e n -V n ), 

(c, + sq 2 ei • e n 



V n ), 



(F.l) 



where V m i n denotes the minimum velocity jump size (due to the cutoff for slow collisions). For V > V m i n , the 
measure on jumps Ve n is 



? m (V,e n ,Q,V)=\ m R 



, 1 M + m 



1m 



V 



ei -e n 



J ni 



1 



ei -e n 



M + m 
1m 



V + V n 



sq 2 



Using this jump measure, we define a Markov process as in Definition IA. 31 
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We write the generator for the Markov process 



L m ^ = V • V Q ^ + / / [V>(Q, V + ^e n ) - V>(Q, V)] f m (V, en, Q, V) dV dfi. 

JS 2 J Vmin 

From the scaling of bath atom velocities and associated law (IF.ll) . E m (|v|) is 0(m~ 1 / 2 ) and by (|F.2[) we have 
that E m (V) is 0(m x l 2 ). We expand the generator for the Markov process in powers of V and find to leading 
order 



L,„r V-V ( ,r(Q.V)-:- C m ( Ii • V V ^(Q, V) + i/ a : V^(Q, V) ) +(Hm> - I 



(F.3) 



where 



L-m — Ait' 



,, , M + m 



M + m 



2m 



i- rco f>2 / i 

A = m- 5 / 2 / / e n , -f m 

JS 2 JVmin l ei ' en l V e l ' e n 

h = m-V 2 / / (e n ® e n ) 1 f / m 



v + v„ 



-s 92 dVdfl, 



1 



M + m 

2m 



K + K 



sq 2 )dVdn, 



where q 2 — Q 2 — i?(e n • ^2). Note that estimating the remainder term follows as in Appendix IB . 1 . 1 1 
In both integrals I\ and I 2 , we substitute 



1 



|ei -e„| 






1/2 e l ' e n 

m ' s- -q 2 , 

ei -e„ 



|ei -e„| 
where x m - m > is the minimum value of x. The drift term, 1%, becomes 

2 



h =m~ 1/2 / / e n (|ei ■ e n \x - m 1/2 V n + m 1/2 s(e 1 • e„) 

JS 2 Va: min V 

8m- 1 / 2 
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» W l ei ' e "l x i dx ^ 

(M + m)3 / Ui-e n Xl ' 



M 3 



s 2 



e n (|ei • e n |a;) / I —a; ) dxdfl 

ei • e„ 



2 Li Gn ( _ml/2K + ml/2s ( e i • e-Ofc) (|ei • enla;) / ( g 1 ^"^ ) dxdn 

/ e„ (-V • e„ + s(e 1 ■ e n )q 2 ) \e ± ■ e n | dO 
VS 2 



+ 0{m 3 ' w ) 



16$i 



+ 0(m 3 / 10 ), 



(F.4) 



(F.5) 



where we have used the fact that 0(m -1 / 2 )-term integrates to zero and have extended the integration region 
with the estimate J* Xmm |ei • e n \xf(x)dx < m x l 2 c m J R f(x) dx = 0(m 3 / 10 ). Note that q 2 = Q 2 — Re n -e 2 depends 
on e„, but the e„-dependence annihilates as we show below in (|F.7|) . 
Substituting x into I 2 , we compute the diffusion coefficient to be 



16 

M^Js 2 
16$ 3 
M 4 



(e n $ e n ) (|ei ■ e n |x) 3 / f ^—^-r X )dxdfl + (^(m 1 / 2 1 

ei • e n | 

l/2\ 



/ (e„®e„)|ei ■ e n \ 3 dfl + 0{m 

J§ 2 



(F.6) 



where we have used the inequality /^"""(lei • e n |a;) 3 /(x) dx < m 3 / 2 ^ J R f(x) dx = 0(vn}l 2 *). 
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Using the following expressions for the spherical integrals 

/ / (e n ® e n )|ei ■ e„| sm<j>d(j)dd 
Jo Jo 



7T 



(e n <8> e n )|ei • e n | 3 sin (j) dcj) d9 = 






77 

2 











71 

2 


2tt 
3 










7T 

6 













/ / (e n <E>e n )|ei • e n |(e 2 • e„) sin^ d<j)d6 = 0, 
Jo Jo 



(F.7) 



we deduce from (IF. 51) and (IF.6|) that 

1 



h = — 



16tt$i 

M3~ 



i (V - s(ei ® ea)Q), 



i 



M4 l I 



Combining these terms with (|F.3|) and (|F.4|) . writing the matrices in terms of tensor products, choosing / 
from (|2.4[) . and taking the limit asm->0, the generator converges in the sense of (|B.2[) to the generator 



£t/; = V-VqV> 



2V2^Ai? 2 



ei ® ei + -(/ - e x ® ei) ) (V - s(ei <g> e 2 )Q) • VyV 



+ BM \3 ei ® Gl + 6 (/ ~ Gl ® ei) ) : V ^^ 
This is the generator for the anisotropic SDE 

dQ = Vdt 

MdV = - 7 (V - s(ei <g) e 2 )Q)d£ + adW = - 7 (V - AQ)dt + adW 



(F, 



where 



7 



2V2ttXR 2 
7$ 



1 
\ 
i 



and 



i 

6 



//3 3 



1/2 



This is the dynamics (J2.12I) discussed in Section 12.31 As mentioned there, we find the inconsistency of (|F.8|) 
when s = with the Langevin dynamics for zero background flow (|2.1[) unacceptable. As an attempt to remove 
the anisotropy observed in (|F.8|) . we add laminar flows in each coordinate direction in the following section. 

F.2. Superimposed baths 

We now consider having three different laminar flows, one for each coordinate direction e^, for % = 1, 2, 3. We 
define the probability distributions 

fi{v) = Zr 1 exp (-^2« 2 J , where Z % = V^9 U i = 1, 2, 3, 

where we allow for different temperatures for each of the bath measures. The initial condition of the bath atoms 
is chosen according to the measure 



d/im(q,v) 



A™™ 1 / 2 



(fi{m l ' 2 {v x - sq 2 ))8(v 2 )S(v 3 ) 



+f 2 (m 1 ? 2 v 2 )S(v 1 )6(v 3 ) + f 3 ('m 1/2 v 3 )8(v 1 )d(v2)) dvdq, 
where 5(x) is the Dirac distribution. 



q,ve 



40 



Upon splitting the integral linearly, the calculations of Appendix IF. II give us the limiting generator 



where 



Lip = V • VqV - M-\jV - u) ■ Vy^ + -M^ 2 a<7 T : AyV 




1/13 
6 W 1 



I #3 





2fl3 
3 W 2 





|6»i 



¥1 






2^2 



^1 






¥1 



¥1 



The corresponding SDE is 



dQ = Vdt, 
MdV = -iVdt ■ 



~idt 



rrfW. 



As an example, let us take 



-kg for i = 1, 2, 3. Then we have for the limiting equation 



dQ = Vdt, 



MdV= - 1 [V--sQ 2 e 1 )dt 



1/2 



rfW, 



where 7 



4Afl J V27r 

3V? 



. This is (|2.15p . which is discussed in Section 



This system satisfies the fluctuation- 



dissipation relation only for inverse temperature j3 — 2/3. That is, the temperature of the large particle is half 
that of the bath. More problematically, there is a factor of 1/2 on s, which means that the large particle only 
feels half the average velocity of the heat bath. This is caused by the superposition of multiple laminar flows. 
We note that there is no choice of 9i,02, and 63 so that the resulting SDE is both isotropic and has response u 
equal to the input background motion. 
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